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ABSTRACT 

We present ultraviolet (UV) observations of six nearby Type la supernovae (SNe la) 
obtained with the Hubble Space Telescope, three of which were also observed in the 
near-IR (NIR) with Wide-Field Camera 3. UV observations with the Swift satellite, as 
well as ground-based optical and near-infrared data provide complementary informa¬ 
tion. The combined data-set covers the wavelength range 0.2-2 /im. By also including 
archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to 
E{B— V) = 1.4 mag. We study the wavelength dependent extinction of each individual 
SN and hnd a diversity of reddening laws when characterised by the total-to-selective 
extinction Ry. In particular, we note that for the two SNe with E{B — V) > 1 mag, 
for which the colour excess is dominated by dust extinction, we find Ry = 1-4 ± 0.1 
and Rv = 2.8 ± 0.1. Adding UV photometry reduces the uncertainty of fitted Ry 
by ~ 50 % allowing us to also measure Ry oi individual low-extinction objects which 
point to a similar diversity, currently not accounted for in the analyses when SNe la 
are used for studying the expansion history of the universe. 

Key words: ISM: dust, extinction - supernovae: general - circumstellar matter 


1 INTRODUCTION 


Studies of the cosmological expansion history using Type la 
supernovae (SNe la) have greatly improved our understand¬ 
ing of the Universe. While pioneering work lead to the dis¬ 
covery of the accelerated expansion from a few dozen SNe la 
i Riess et al.||1998 Perlmutter et al!]|1999 l the present sam¬ 
ples of several hundred SNe la out to 2 ~ 1.5 show to high- 
precision that Einstein’s Cosmological Constant provides an 


excellent fit for this phenomenon (e.g. Amanullah et al. 2010 
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20141. 

idings have been possible due to the great ac- 
e la as distance indicators when observed in 
broadband filters over several weeks, starting 
urve peak (see Goobar & Leibundgut|20111 for 
w). The distance to each individual SN is ob- 
ing a SN la template to the photometric data 

(see e.g. Riess, Press & Kirshner|1996 Nugent, Kim & Perl- 

mutter||20021 

Jha, Riess & Kirshner||2007t |Guy et al.||20071 

Conley et al.| 

2008 Burns et al. 12011 and references therein) 

in order to obtain the brightness and colour at maximum to¬ 
gether with the light curve-shape. The fitted properties are 


then typically combined using empirical linear relations to 
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form a distance-dependent quantity (Phillips 1993 Tripp 


1998). 


ing the relation between SN colours and spectroscopic prop- 


The colour correction is based on the observation that 
red SNe are fainter than bluer ones, which could originate 
from a combination of an intrinsic SN la colour law, and 
extinction by dnst in the host galaxy. The colour-brightness 
relation is commonly determined by minimising the scatter 
aronnd the best fit cosmological model, but there is no a 
priori reason why the two effects should follow the same 
colour-luminosity relation, nor can it be expected that host 
galaxy dust should have the same properties in different SN 
environments. In fact, even in the Milky Way, a range of 
extinction properties have been observed for different lines- 
of-sight. 

When SNe la are used to measure distances, these un¬ 
certainties are handled by introducing a systematic error 
on the colour-luminosity relation, which makes an impor¬ 
tant contribution to the error budget when propagated to 
the derived cosmological parameters. Addressing the colour¬ 
brightness relation, and in particular breaking the degener¬ 
acy between the intrinsic colour component and dust ex¬ 
tinction, is important for future SN la surveys to improve 
beyond the current cosmological constraints. 

The properties and wavelength-dependent extinction of 
dust in the Milky Way has been carefully studied and is 
commonly characterised by the total, Ay, to selective, E{B— 
V) — Ab — Av, extinction ratio as Rv = Av/{Ab — Ay). 
Lower values of Ry correspond to steeper extinction laws, 
since, for a given total extinction. Ay, these imply a larger 
reddening, E{B — V). 

The reddening of SNe la can be studied by comparing 
observed colours between reddened and similar un-reddened 


objects. Several such studies of individual SNe la (e.g. Elias- 


Rosa et al.|2006||2008[[Krisciunas et al.|2007||Folatelli et’ah 


2010[ Amanullah et al. 2014[ from hereon A14) in nearby 
galaxies, point to significantly lower values of Ry than has 
been observed in the Milky Way. Since the reddening of 
several of these SNe la is significantly higher than what can 
be expected from the Milky Way extinction along the line of 
sight, the observed colour excesses are likely dominated by 
extinction in the SN host galaxies. For low reddening SNe 
this may not be the case, but it is by nature difficult to 
measure Ry to high precision for individual SNe with low 
E{B— V). However, the global reddening law for a sample of 
SNe can be measured, and several studies have obtained low 
Ry values (e.g. Astier et al. [2006 Nobili & Goobar 2008 


Burns et al.|2014 from hereon B14) using this approach. 


erties. Wang et al. (2009) found Ry — 1.6 for SNe la with 


’’high” photospheric velocities (HV, > 11,800 km/s) within 
5 days of maximum, while they obtained Ry = 2.4 for ob¬ 
jects with ’’normal” velocity (NV, < 11,800 km/s). The ex¬ 
pansion velocities were quantified by measuring the veloc- 


ity of the Si ll 6355 A absorption feature. Foley & Kasen 


(2011) confirmed their findings but also conclude that when 


SNe with E{B — V) > 0.35 mag are omitted, they obtain 
Ry ~ 2.5 for the two separate subsamples, which also have 
different intrinsic colours. In other words, these findings also 
suggest that low Ry values are primarily associated with 
high-extinction, and these could dominate when a single 
global reddening law is derived from a sample of SNe la. 

One natural explanation for low Ry values arises if 
SNe la are surrounded by circumstellar (CS) dust (|Wang| 


2005 Goobar 2008). In this scenario, photons can scatter 


back into the line-of-sight which reduces the total extinc¬ 
tion. Photons that scatter on the CS dust will also arrive 
later than photons that do not interact which will give rise to 
time-dependent reddening ( Wang|2005 Amanullah & Goo- 


barpOlf Brown et al.pOiS ). An observation will be a su¬ 

perposition of the SN spectrum at the given epoch and SN 
light from earlier epochs, but the scattering cross-section 
is wavelength-dependent so the fraction of delayed photons 
will also vary with wavelength. Further, depending on the 
geometry, scattered photons could give rise to a plateau for 
the late-time tail of the light curve and indications of such 
tails have been observed in B — V ( Forster et al.p013 |. 

CS dust could also be heated by the SN, where the tem¬ 
perature will depend on the dust properties and the distance 
from the explosion. Heated dust has been detected in mid 


to far-infrared for a subset of peculiar SNe la ( 

Fox et al. 

2011||Taddia et al.|2012||Fox & Filippenko|2013t 

Silverman 

et al.|2013 

1 , but Herschel and Spitzer observations ( 

Johans- 

son, Amanullah & Goobar 2013 Johansson et al.||2014l of 


a handful of normal SNe la show no signs of heated CS 
dust (see also Maeda, Nozawa fc MotoharapOfi for limits 
from Near-IR observations). For two of these SNe, 2006X 
and 20071e, together with SNe 1999cl and the SN 2002ic- 
like PTFllkx, time-varying Na l D absorption features have 
been observed in their spectra on time-scales comparable 
to the SN life-time which is consistent with the presence 
for CS material JPatat et al.||2007 jBlondin et al/]|2009 Si 


mon et ar]|2009 [Pilday et al.p012i. Absorption from Na 

has been used as a proxy for dust extinction (see e.g. |Poz-| 
nanski, Prochaska fc Bloom|2012'|, although high-resolution 


that globally derived reddening laws will be biased if the 
intrinsic SN la colour dispersions are not accounted for cor- 


Chotard et al. (2011) and Scolnic et al. (2014) noted spectra are required (Poznanski et al. 2011) and its valid- 


rectly, and Chotard et al.l (20111 find Ry = 2.8 ± 0.3, con¬ 


sistent with the Milky Way average, for their sample of 
76 SNe la of which 73 had E{B — V) < 0.3 mag. B14 also 
obtain Ry ~ 3 when they constrain their analysis to only 
include low-extinction SNe la, while highly reddened ob¬ 
jects appear to prefer lower values {Ry ~ 1.7) which is in 
agreement with previous studies ( Folatelli et al.|2010 l. This 
suggests that a single, global, reddening law derived from a 
set of SNe la could depend on the colour distribution of the 
sample. 

The degeneracy between intrinsic SNe la colour varia¬ 
tions and dust extinction, can also be approached by study- 


ity for extinction along lines-of-sight of SNe la has been 
challenged ( Phillips et al.p013 l. Time-varying Na l D can 
be explained by photo-ionisation of neutral sodium (e.g. 
[Borkowski, Blondin fc Reynolds||2009| ) and subsequent re¬ 
combination by the SN. In a recent study by Soker (2014) 


this explanation has been questioned and an alternative 
model of photon-induced desorption of sodium from dust in 
planetary nebulae remnants was proposed. It has also been 
suggested that time-variations may originate from patchy in¬ 
terstellar material along the line-of-sight ( Patat et al.|2010 1. 
More recently, Graham et al. (2015b I claimed a detection of 


a time-varying K l absorption feature in 2014J, while the 
corresponding Na l D feature remained unchanged, which 
they argue is consistent with the presence of CS material at 
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a radius of rdust ~ 10^® cm due to the different ionisation 
cross-sections. 

An extensive study with multi-epoch high-resolution 
spectra of 14 SNe ( Sternberg et al.|2014 | have failed to re¬ 
veal further examples of time varying Na l D. Nevertheless, 
a statistically significant preponderance of Na l D features 
which are blue shifted with respect to the local velocity have 


been observed in high- and mid-resolution spectra (Stern- 
berg et al.|2011 Maguire et al. respectively), suggest¬ 


ing that there is outflowing material from the SN progenitor 
system. 

In this work we seek to expand our knowledge of ex- 
tragalactic extinction in the line-of-sight and SN la redden¬ 
ing by making the widest wavelength study of a sample of 
SNe la to date. We present UV-to-NIR light curves for six 
SNe la based on data mainly obtained from The Nordic 
Optical Telescope (NOT) and The Hubble Space Telescope 
(HST). The observations and data reduction is presented in 
§[^ while the light curves and colours are presented in §|^ 
The method for deriving individual reddening laws for each 
SN is described in §|^ and after adding SN 2014J from A14 
to the sample we present the results in § |4.3| We search for 
CS dust in §[^ and all the results are discussed in §|^ We 
conclude and summarise in §[^ 

All seven SN^are listed in Table and shown together 
with their host galaxies in Figure They are also briefly 
summarised in Appendix [X] 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Hubble Space Telescope 

All previously unpublished HST observations discussed in 
this work are listed in Tabled Each SN was observed at two 
different epochs (four for SN 2012et) with the Wide-Field 
Camera 3 (WFC3) during the HST Cycle 19 under pro¬ 
gramme GO-12582 (PI: Goobar). All SNe were imaged with 
the WFC3/UVIS channel through the passbands F225W, 
F275W and F336W. In addition to this, NIR imaging in the 
HST passbands F125W and F160W was obtained for the 
SNe 2012cg, 2012cu and 2012et with the WFC3/IR channel. 
Examples of both the WFC3/UVIS and WFC3/IR observa¬ 
tions are shown in Figure 

WFC3/UVIS consists of two 4096 x 2051 e2V CCD 
detectors with a plate scale of 0.04”/pixel, while the Tele¬ 
dyne HgCdTe infrared detector used in the WFC3/IR chan¬ 
nel has a pixel scale of 0.13”/pixel. For each observation, 
only the sub-arrays (apertures), listed in Tableof the de¬ 
tectors were read-out. All WFC3 data were reduced using 
the standard STScI reduction pipeline and calibrated using 
CALWF3. 

The WFC3/UVIS CCDs, like all HST CCDs, are 
plagued with charge transfer inefficiencies (CTI), i.e. degra¬ 
dations of the detector performance over time due to damage 
in the silicon lattice from cosmic rays. The CTI can be re¬ 
verse corrected at the pixel level using wf c3uv_cterevers^ 


^ All data and figures presented in this work are available at 
http://snova.fysik.su.se/dust/ 

^ http://WWW.Stscl.edu/hst/wfc3/tools/cte_tools 


Civil date 

MJD 

Aperture 

SN 

2012-04-09 

56026.4 

UVIS1-2K2A-SUB 

2012bm 

2012-04-09 

56026.5 

STIS 52X0.2 G230LB 

2012bm 

2012-04-13 

56030.6 

UVIS1-2K2A-SUB 

2012bm 

2012-04-13 

56030.7 

STIS 52X0.2 G230LB 

2012bm 

2012-04-16 

56033.1 

STIS 52X0.2 G230LB 

2012bl 

2012-04-16 

56038.2 

UVIS1-2K2A-SUB 

2012bl 

2012-04-20 

56037.4 

STIS 52X0.2 G230LB 

2012bl 

2012-04-21 

56038.2 

UVIS1-2K2A-SUB 

2012bl 

2012-06-04 

56082.4 

UVIS2-M1K1C-SUB 

2012cg 

2012-06-04 

56082.4 

IRSUB64 

2012cg 

2012-06-04 

56082.5 

STIS 52X0.2 G230LB 

2012cg 

2012-06-04 

56082.5 

STIS 52X0.1 G430L 

2012cg 

2012-06-18 

56096.4 

STIS 52X0.2 G230LB 

2012cg 

2012-06-18 

56096.4 

STIS 52X0.1 G430L 

2012cg 

2012-06-19 

56097.5 

UVIS2-M1K1C-SUB 

2012cg 

2012-06-19 

56097.5 

IRSUB64 

2012cg 

2012-06-04 

56082.6 

UVIS1-2K2A-SUB 

2012cp 

2012-06-04 

56082.7 

STIS 52X0.2 G230LB 

2012cp 

2012-06-16 

56094.6 

UVIS1-2K2A-SUB 

2012cp 

2012-06-17 

56095.6 

STIS 52X0.2 G230LB 

2012cp 

2012-07-02 

56110.8 

UVIS2-M1K1C-SUB 

2012CU 

2012-07-02 

56110.8 

IRSUB512 

2012CU 

2012-07-07 

56115.8 

UVIS2-M1K1C-SUB 

2012CU 

2012-07-07 

56115.8 

IRSUB512 

2012CU 

2012-10-01 

56201.5 

UVIS2-M1K1C-SUB 

2012et 

2012-10-01 

56201.5 

IR 

2012et 

2012-10-01 

56201.8 

STIS 52X0.2 G230LB 

2012et 

2012-10-01 

56201.8 

STIS 52X0.2 G430L 

2012et 

2012-10-05 

56205.6 

UVIS2-M1K1C-SUB 

2012et 

2012-10-05 

56205.6 

IR 

2012et 

2012-10-09 

56209.8 

UVIS2-M1K1C-SUB 

2012et 

2012-10-09 

56209.8 

IR 

2012et 

2012-10-13 

56213.9 

UVIS2-M1K1C-SUB 

2012et 

2012-10-13 

56213.9 

IR 

2012et 


Table 2. Wide-Field Camera-3 and STIS observations together 
with the used apertures, where UVISl and UVIS2 are the two 
chips of the WFC3/UVIS channel. The observations with the 
WFC3/IR and WFC3/UVIS channels were obtained through the 
F125W, F160W and F225W, F275W, F336W broadband pass- 
bands respectively. 


for all frames except for images obtained with the UVIS2- 
MIKIC-SUB sub-array. For this aperture we instead used 
the recipe suggested by Bourque & Anderson (20141. The 
individual flat-fielded WFC3/UVIS and WFC3/IR images 
were then resampled and combined using astrodrizzle 
( Fruchter fc et al.|2010 |. 

Aperture photometry was carried out on the combined 
images using an aperture radius of 0.4”, although in some 
cases a smaller aperture of radius 0.16” had to be used 
for the WFC3/UVIS data. This was the case for the low 
signal-to-noise measurements of SN 2012bm, observations of 
SN 2012cg that were plagued by nearby coincidence cosmic 
rays, and all epochs of SN 2012et, where a nearby object 
contaminated the larger aperture (see the mid-right panel 
of Figure]^. The smaller radius measurements had to be 
aperture corrected to the 0.4" radius for which STScI pro- 
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Figure 1. The seven SNe analysed in this work. The patches for SNe 2012bm, 2012cg, 2012cp, 2012cu and 2012et were obtained 
from BVR images from NOT, while the images for SNe 2012bl and 2014J are obtained from the Swope telescope and HST/WFC3 
(Program DD-13621; PI Goobar) images respectively. 


SN 

“SN 

(2000) 

(SSN 

(2000) 

r± 

(kpc) 

Host galaxy 

“h 

(2000) 

<5h 

(2000) 

(km/s) 

Distance 

(Mpc) 

MW Av 
(mag) 

2012bl 

20:23:55.28 

-48:21:17.3 

15.6(0.2) 

ESC 234-19 

20:23:51.0 

-48:21:32 

5608 

71(1)“ 

0.098 

2012bm 

13:05:45.66 

-1-46:27:52.9 

7.5(0.3) 

UGC 8189 

13:05:46.6 

+46:27:42 

7436 

103(8)*' 

0.033 

2012cg 

12:27:12.83 

-1-09:25:13.2 

1.3(0.2) 

NGC 4424 

12:27:11.6 

+09:25:14 

437 

15(2)“ 

0.057 

2012cp 

13:47:01.79 

-1-33:53:35.0 

1.9(0.1) 

UGC 8713 

13:47:01.2 

+33:53:37 

4956 

57(3)** 

0.058 

2012CU 

12:53:29.35 

-1-02:09:39.0 

5.4(1.2) 

NGC 4772 

12:53:29.1 

+02:10:06 

1040 

41(9)“ 

0.074 

2012et 

23:42:38.82 

-1-27:05:31.5 

2.7(0.2) 

MCG -1-04-55-47 

23:42:38.4 

+27:05:31 

7483 

105(7)“ 

0.221 

2014J 

09:55:42.11 

-1-69:40:25.9 

1.0(0.3) 

M82 

09:55:52.7 

+69:40:46 

203 

3.5(0.3)f 

0.189f 


Parnovsky &: Parnowski 


Springob et al. l[200ij I 


I 2010 1; ^ Calculated based on the redshift with Ho = 73 it 5 km/s/Mpc; Cortes, Kenney & Hardy I 20081 


Itilly et al. (2009 i; ^ Dalcanton et al. (20091 


Table 1. SN coordinates are quoted from the discovery telegrams. All host galaxy data were obtained from the Nasa Extragalacic 
Database (NED), unless otherwise specified, where v^^ is the measured recession velocity, and the Milky Way extinctions are from the 
[Schlafly &: Finkbeiner| ( |2011| | calibration of the |Schlegel, Finkbeiner &: PavT^ | |1998| l infrared-based dust maps. The projected distances 
from the host galaxy nuclei, ry, were calculated based on the host galaxy distances and the SN offsets specified in Appendix [X] 


vides zeropoint^for WFC3. Since there were no bright stars 
present in the fields we used the general enclosed energy 
tables provided by STScI for this purpose. To account for 
the time and wavelength dependence of the WFC3/UVIS 
aperture corrections we adopt an additional uncertainty of 
0.05 mag for all measurements using the 4 pixel radius, 
which was derived by comparing the enclosed energy tables 


® We used the zeropoints from March 6th, 2012. 


with photometry from different aperture radii of the high 
signal-to-noise data. 

Only the SN 2012et NIR measurements suffered from 
host galaxy contamination. The host contribution was es¬ 
timated by placing apertures of the same radius used for 
the SN photometry along the isophot of the galaxy that in¬ 
tersects with the position of SN 2012et as illustrated in the 
lower right panel of Figure The background was then es¬ 
timated as the median of these measurements, and the root- 
mean-square of the measurements was added in quadrature 
to the photometric uncertainty. All the WFC3 photometry 
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SN2012bl 


SN2012bm 


SN2012cp 

F275W 


F275W 


F275W 


SN2012cg 


F275W 


- • 


F160W 


SN2012CU 

F275W 




F160W 




Figure 2. HST/WFC3 patches of the six SNe observed during 
Cycle 19. The patch sizes are 10" X 10" and in the upper row the 
first epochs obtained through the F275W filter is shown for each of 
the three SNe 2012bl, 2012bm and 2012cp. In the two lower rows, 
the corresponding patches are shown for the SNe 2012cg, 2012cu 
and 2012et, together with the NIR observations in the F160W 
filter that also were obtained for these objects. The marked (red) 
region in the lower right panel show the isophot from which the 
host galaxy background of SN 2012et was estimated as explained 
in the text. 


is presented in Tableand shown in Figure]^ Here we also 
present updated photometry of SN 2014J using the method 
described here which is consistent with the results presented 
in A14. 

For most HST visits we also obtained long-slit spec¬ 
troscopy with the Space Telescope Imaging Spectrograph 
(STIS) where the 1024 x 1024 pixel SITe CCD detector 
was used with the G230LB grating covering the wave¬ 
length range 1600-3100 A, and for two SNe (2012cg, 2012et) 
G430L, 2900-5700 A. The data were reduced using the 
calstis pipeline which is part of the STSDAS package. The 
pipeline was ran for all the spectra up to the point where 
the calibrated 2D spectra were created. Only the spectra of 
SN 2012cg and the red G430L spectrum of SN 2012et con¬ 
tained any significant signal and were extracted using the 
calstis xld routine with a 4 pixel aperture. 


2.2 Swift/UVOT 


UV photometry was also obtained with the Ultra- 


Violet/Optical Telescope {UVOT, Roming et al. 20051 on 
the Swift spacecraft ( Gehrels et al.|2()04 1 for the SNe 2012bl, 
2012cg, 2012cp and 2012cu in the uvwl, uvw2 and uvm2 fil¬ 
ters. However, the uvwl and uvw2 filters are not well suited 
for extinction studies due to the significant “red-tails” of 


these filters. Milne et al. (20101 used the UV spectra of 


SN 1992A (Kirshner et al. 1993 \ and estimated that 52 % 


and 44 % of the light in uvwl and uvw2 resp ectively, orig- 
inate from wavelengths redder than 3000 A. Brown et al. 


(20151 show how the observed photons and the correspond¬ 
ing effective wavelengths shift dramatically as reddening in¬ 
creases. 

The uvm2 filter, on the other hand, is much better con¬ 
strained (1 % of the light comes from > 3000 A) and is used 
for the extinction studies in this work. The uvm2 magni¬ 
tudes were measured using the pipeline from the Swift Op¬ 


tical/Ultraviolet Supernova Archive (SOUSA, Brown et al. 


20141, including correction for the time-dependent sensitiv¬ 
ity and revised zeropoints from [Breeveld et al. (20111. The 
result is presented in Table and shown in Figure [3 Only 
the SNe 2012bl and 2012cg showed a significant SN flux 
in the uvm2 band after subtraction of the underlying host 
galaxy flux. The flux measured for SN 2012cp was consistent 
with the host galaxy brightness and no signal was measured 
for SN2012CU. 


2.3 Ground based observations 


The ground based spectroscopic observations are listed in 
Table 1^ while the photometric observations are summarised 
together with the measured magnitudes in Table and 
shown in Figure All SNe except for SN 2012bl were ob¬ 
served with the 2.56 metre Nordic Optical Telescope (NOT) 
under programmes 45-009 and 46-018 (PI: Amanullah). 
Imaging and spectroscopy was obtained with the 6.4' x 6.4' 
Andalucia Faint Object Spectrograph and Camera (AL- 
FOSC) using the filter set U (#7fl B (#74), U (#75), R 
(#76), i (#12) and the R — 360 grism (#4)[^ The data 
were reduced using standard IRAF routines and the QUBA 
pipeline ( [Valenti et al. 20111. The photometry was mea¬ 
sured by first fitting the point-spread-function (PSF) using 


DAOPHOT (Stetson 19871 to stars in the fields. The cali¬ 


brated magnitudes were obtained either from Landolt fields 
(Landolt 19921 observed during the night, or, when these 


were not available, by comparing the fluxes to stars in the 
fields that were then calibrated against Landolt fields for the 
photometric nights. 

With the NOT we also obtained high-resolution spec¬ 
troscopy of SN 2012cg using the Fibre-fed Echelle Spectro¬ 


graph (FIES, Telting et al. 20141 in its high-resolution mode, 
R = 67000. Simultaneous wavelength reference (Thorium- 
Argon) spectra were obtained and the data were reduced 
using the software FIEStool which is provided by the obser¬ 
vatory. 

NIR observations were carried out for SNe 2012bm, 
2012cg and 2012cp with the 4' x 4' NOTCam instrument 
in the J, H and Ks bands under the programmes 45-007 
(PI: Kankare) and 46-020 (PI: Mattila). We used the wide 
field imaging option of the 1024 x 1024 pixel HgCdTe NOT¬ 
Cam detector with a plate scale of O'.'234/pixel. The obser¬ 
vations were carried out using either 5 or 9-point dithering 
patterns. For the SNe with extended host galaxies, beam¬ 
switching was used to guarantee a successful sky subtraction. 
We used the NOTCam Quick-Look reduction package based 
on IRAF. Bad pixels (which includes two dead columns) 
were masked and we used master differential skyfiats which 


4 NOT filter ID 
® NOT grism ID 
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Figure 3. Observed light curves for the SN sample. Our analysis is based on optical UBVRI and NIR JHK photometry from ground 
based observatories (cirles), UV (F225W,F275W,F336W) and NIR (F125W,F160W) photometry from the HST (squares). For the SNe 
2012bl and 2012cg we also use Swift/ UVOT photometry (triangles). The data are shown in their instrumental system but have been 
shifted using the distance moduli based on the distances in Tablej^and offset for clarity. The thick, solid, lines show the SNooPy F-models 
used to calculate colours as described in the text, while the dashed lines show synthesised light curves of SN 2011fe to guide the eye. 
When a SN has observations in a similar filter from multiple instruments only one of the light curves is shown. The epoch is the observed 
date with respect to the fitted day of maximum in the R-band from Table 


were obtain from bright and faint skyflats. Further, the im¬ 
ages were corrected for geometric distortion^ before the in¬ 
dividual images were aligned and coadded. 

SN 2012cg was also observed in JHKs with the 4.3' x 4.3' 
CAIN instrument on the 1.52 metre Carlos Sanchez Tele¬ 
scope at Observatorio del Teide on Tenerife. CAIN III is a 
256 X 256 pixel HgCdTe NIR detector, where we used the 
wide field option with a plate scale of 1.0"/pixel. The data 
were reduced using a dedicated IRAF package provided by 
J. Pullido and A. Barrena. 

The two SNe 2012bm and 2012cp were also observed 
by the Palomar Transient Factory with the Palomar 48- 


inch telescope (Ran et al. 20091 in the Mould R. We use 


® Magnus Galfalk, private communication. 


their photometry which was obtained by first determining 
the PSF prior to subtractions, and then carrying out PSF 
photometry on the subtracted frames | Firth et al.| 2015|. 
The calibration of this data has been described in lOfek et al.l 


( 20121 . 


Photometric follow-up of SN 2012bl was carried out by 
the Carnegie Supernova Project using the Swope telescope 
for the optical and the du Pont telescope for the NIR. The 
reduction and photometry have been described in detail in 


Hamuy et al. (20061 and Contreras et al. (20101. PSF pho¬ 


tometry was performed with respect to a local sequence of 
standard stars calibrated to the [Landolt] ( |1992[ ) and |Smith| 
et al. (20021 standard fields for the optical observations and 


the Persson et al. (19981 for the NIR. 


Further, we here also present spectroscopic observations 
of SNe 201 Ife, 2012cg and 2014J obtained using the ISIS 
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Diversity in extinction laws for SNe la 7 


MJD 

Phase 


A 

^MW 

Kx 


V 

^MW 

Kv 

o 

1 

(days) 

(days) 

Filter 

(mag) 

(mag) 

(mag) 

Match 

(mag) 

(mag) 

(mag) 

(mag) 

SN2012cg 

56082.4 

0.2 

WFC3 F336W 

11.74(0.06) 

0.09 

0.01 

M 

12.04(0.04) 

0.06 

0.00 

0.61 

56097.5 

14.0 

WFC3 F336W 

13.75(0.01) 

0.09 

0.02 

M 

12.53(0.07) 

0.06 

0.00 

-0.90 

56082.4 

0.2 

WFC3 F125W 

12.37(0.01) 

0.01 

0.05 

M 

12.04(0.04) 

0.06 

0.00 

-0.55 

56097.5 

14.1 

WFC3 F125W 

13.88(0.03) 

0.02 

-0.07 

M 

12.53(0.07) 

0.06 

0.00 

-1.91 

56082.4 

0.2 

WFC3 F225W 

16.38(0.05) 

0.13 

0.02 

M 

12.04(0.04) 

0.06 

0.00 

-2.55 

56097.5 

14.0 

WFC3 F225W 

17.60(0.08) 

0.13 

0.02 

M 

12.53(0.07) 

0.06 

0.00 

-3.68 

56082.4 

0.2 

WFC3 F160W 

12.67(0.02) 

0.01 

0.01 

M 

12.04(0.04) 

0.06 

0.00 

-0.88 

56097.5 

14.1 

WFC3 F160W 

13.24(0.03) 

0.01 

0.55 

M 

12.53(0.07) 

0.06 

0.00 

-0.69 


Table 3. The photometry of all SNe. All magnitudes are in the natural Vega system. The rest-frame magnitude and colour can be 
obtained from Columns 4—6 using Equation 0 - Columns 5 and 9 are the Galactic extinctions and Columns 6 and 10 are Kx and Ky 
corrections for the two bands, respectively. All corrections were calculated after the SN 2011fe template had been reddened with the best 
fitted F99 law, shown in Table[^ for each SN. The V magnitude is only shown for data points used in the colour analysis, i.e. with phases 
between —10 and +35 days. Column 2 show the effective light curve-width-corrected phase, while Column 7 specifies whether the V 
magnitude was measured for the same epoch (D) or if it was calculated using the SNooPy model (M). The corresponding intrinsic colour 
for the SN 2011fe template is shown when available in Column 11. The Swift magnitudes of SN 2014J are from |Brown et ahjpOlSl . (This 
table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form 
and content.) 


Civil date 

MJD 

Epoch 

(days) 

Instrument 

A 

(A) 

Exp. 

(s) 

SN 

2012-05-26 

56073.8 

+260 

WHT/ISIS 

3500-9500 

600 

20 life 

2012-03-30 

56016.3 

-2 

DuPont/WFCCD/WF4K-l 

3600-9200 

700 

2012bl 

2012-03-31 

56017.4 

-1 

DuPont/WFCCD/WF4K-l 

3600-9200 

700 

2012bl 

2012-05-01 

56048.3 

+29 

DuPont/WFCCD/WF4K-l 

3600-9200 

900 

2012bl 

2012-03-28 

56014.0 

-4 

ASIAGO/AFOSC 

3500-8200 

2700 

2012bm 

2012-04-09 

56027.0 

+9 

NOT/ALFOSG 

3200-9000 

2400 

2012bm 

2012-04-16 

56034.1 

+16 

NOT/ALFOSG 

3200-9000 

2400 

2012bm 


2012-05-24 

56071.9 

-9 

NOT/ALFOSG 

3700-7200 

5477 

2012cg 

2012-05-26 

56073.9 

-7 

NOT/FIES 

3640-7360 

1800 

2012cg 

2012-05-28 

56075.9 

-5 

NOT/ALFOSG 

3700-7200 

6588 

2012cg 

2012-06-02 

56080.9 

+1 

NOT/FIES 

3640-7360 

900 

2012cg 

2012-06-04 

56082.9 

+3 

NOT/ALFOSG 

3200-9000 

300 

2012cg 

2012-06-14 

56092.9 

+13 

NOT/FIES 

3640-7360 

1200 

2012cg 

2013-03-09 

56360.4 

+280 

Keck/LRIS 

3075-10300 

600 

2012cg 


2012-05-26 

56074.0 

-7 

NOT/ALFOSG 

3200-9000 

2700 

2012cp 

2012-06-04 

56083.0 

+2 

NOT/ALFOSG 

3200-9000 

1800 

2012cp 

2012-07-05 

56113.9 

+9 

NOT/ALFOSG 

3200-9000 

1200 

2012CU 

2012-09-13 

56183.8 

-7 

ASIAGO/AFOSG 

3500-8200 

2700 

2012et 

2012-09-22 

56193.1 

+2 

NOT/ALFOSG 

3200-9000 

600 

2012et 

2014-11-09 

56970.0 

+282 

APO/DIS 

3300-9800 

750 

2014J 


Table 4. Ground based spectroscopic observations. The Asiago spectra were obtained through the Asiago Transient Classification 
Program ([Tomasella et al.|2014[|. The epochs are shown with respect to the date of B-maximum shown in Table 


instrument at the William Herschel Telescope, the LRIS in¬ 
strument at Keck and the DIS instrument at the Apache 
Point Observatory. These data were reduced in a standard 
manner. 


3 SN SPECTRA AND COLOURS 


The spectra for the seven SNe around maximum are shown 
in Figure where we also added the classification spec¬ 
trum of SNe 2012bm and 2012et (Tomasella et al. 20141 
to the set of spectroscopic observations already described 
above. For SN 2012cg, both the HST/STIS G250LB and 
G430L spectra, and the ALFOSC spectrum, that was ob- 
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8 R. Amanullah et al. 



The colours were also K corrected (e.g. Nugent, Kim 
& Perlmutter 2002 1 and S corrected (e.g. Suntzeff 2000 


Restframe wavelength, A [A] 

Figure 4. Spectra near maximum, when available, for the SNe in 
our sample. The grey spectra are for SN 2011fe at similar epochs 
( [Pereira et al.|[2013| [Mazzali et al.||2014| | reddened by the ex¬ 
tinction laws that we derive in §|4| For SNe 2012cg and 2014J 
we have combined the optical spectra with STIS spectra (the 
SN 2014J spectrum was taken from [^ley et aL||2014^ . Spectro¬ 
scopic features discussed in §[^ §|6.1|and §|6.4|have been marked. 


tained the same day, are shown. We used SNID (Blondin & 
Tonry|2007 1 to type the spectra and all SNe could be classi¬ 
fied as ’’normal”, although SN2012et also provided accept¬ 
able matches to SN 1999aa. For each spectrum in Figure 
we also show the spectrum of the normal SN 201 Ife for the 
matching epoch ( Mazzali et al.|2014 l for comparison. 

All photometric data are shown in Figure and 
for each measurement in each filter, X, where X € 
{wum.8, F225W,..., As}, we calculated the colour X —V. 
When X and F-band measurements are not available for the 
same dates, and the colour could not be obtained directly, 
a F-band model was used to calculate the colours. When 
applicable, we used a smoothed spline model but for the 
SNe 2012bm and 2012cp the F-band data were too sparsely 
sampled to allow reliable spline fits. For these objects we 
instead fitted the SN la template from [Hsiao et al.j | |2007| 
H07 from hereon) to the observed data. Both the spline and 
template-based models were obtained using the SNooPy light 
curve fitter ( [Burns et al. TMTt and are shown as thick lines 
in Figure [^ Further, for all model-based colours we add 
0.05 mag in quadrature to uncertainties in order to account 
for the inaccuracy of the models. 

All colours were corrected for Galactic extinction, 
using the extinction law by [Cardelli, Clayton fc Mathisj 
(1989 CCM from hereon) with Rv ~ 3.1 and the 


values shown in the last column of Table [T] The values 
are from the Schlafiy & Finkbeiner ( 2011[ ) recalibration of 
the Schlegel, Finkbeiner & Davis 119981 dust maps, except 
for SN 2014J where we use a value obtained from studying 
neighbouring regions (Dalcanton et al.|2009|. 


Stritzinger et al.|2002 Krisciunas et al.|2003 1 to a common 

rest-frame filter system that will be used for the remaining 
analysis throughout the paper. The analysis was carried out 
using the WFC3/UVIS and Swift/ UVOT S\tem for the UV, 
the ALFOSC filter set for the optical and the NOTCam fil¬ 
ters for the NIR, since the bulk of the data were obtained 
using these filters. The combined K and S corrections {Kx 
from now on) are calculated synthetically using the SED of 
SN 20life as described in § |4.1| and the filter transmissions 
provided by the different observatories. 

To summarise, all colours were obtained as 

X-V = {mx-AT'-Kx)-{Mv-A^^-Kv), (1) 

where mx are the measured instrumental magnitudes in fil¬ 
ter X and A4 y is the F magnitude for the same date either 
measured directly, or obtained from a spline or template 
model. All values are presented in Table [^ where we have 
also added the Swift photometry of SN2014J from [Brownj 
et al. (20151 to the measurements described in §|^ As will 


be further discussed in §|^ not only the observed colours, 
but also, and Kx depend on the reddening of each SN 
and were therefore obtained iteratively during the extinction 
law fitting. The values shown in Table [^ were obtained after 
the SN20life SED had been reddened with the jFitzpatrickj 
( 1999[ hereafter F99) law using the parameters shown in 
Table m for each SN. 

A subset of the colours are also shown in Figurej^where 
the colours obtained using measurements and the F-band 
model are shown with filled and open symbols respectively. 

SNe la are intrinsically bright in the rest-frame B-band 
and this filter is commonly used to measure SN la distances. 
SN la light curves are typically quantified by a few parame¬ 
ters such as the brightness and time of maximum, for 

e.g. the B-band, and the light curve shape. The light curve 
shape can be quantified using different methods, where per¬ 
haps the most straight-forward is to measure the brightness 
decline for the first 15 days past Phillips (19931 showed 


that fast declining SNe la are typically intrinsically fainter 
than slow declining objects, a discovery that lay the foun¬ 
dation for using SNe la as distance indicators for precision 
cosmology. An alternative approach to quantify the light 
curve shape is to introduce a stretch parameter ( [Perlmutter[ 
et al. 19971, sb, which can be defined as the value needed 


to match the time-evolution of an observed B light curve to 
a standard SN la template. 

The time-evolution of SN la A — F colours also depends 


on light curve shape. For example, Nobili & Goobar (2008 


from hereon N08) present colour-light curve shape relations 
in the optical while B14 show that the SNe la colours can 
be standardised over a wide range of decline rates by in¬ 
troducing colour-stretch, sbv = tsy/SO days, where tBv 
is the time between and the maximum of the B — F 
colour. However, this approach requires B — F coverage up 
to t ~ t’f/’^^ 40 in order to accurately determine sbv- 

We here standardise the measured colours using s b ■ All 
colours will be studied and compared as a function of phase, 
p, defined as p = (t — t’B^^)/sB, where t are the observing 
dates. The phases for the colours shown in Figure have 
been determined by using the values of and sb from 
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Diversity in extinction laws for SNe la 9 


▲ 

SN2011fe 

T 

SN2012bm 

■ SN2012CP 

► 

SN2012et 

★ 

SN2012bl 

• 

SN2012cg 

◄ SN2012CU 

♦ 

SN2014J 



Figure 5. A selection of measured colours from UV to NIR together with literature SNe 2011by and 2011fe. The phases, p = (f— 
were obtained using the values of and sb from Table [s] The colours from the H07 spectral template (dashed red) and |Nobili &| 

[Goobar] l |2008[ l (NOS, solid blue with dispersion region) are mso shown together with the derived colour model of SN2011fe (solid grey) 
with the adopted dispersions (grey region). Filled symbols show colours where a F-band measurement was obtained for the same date, 
while a F-band model was used to obtain the measurements that are shown with open symbols. Errorbars are only plotted for data 
points where the uncertainty exceeds 0.1 mag. See text for details. 
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10 R. Amanullah et al. 


Table which are fitted simultaneously with the reddening 
law as described in §|^ 


4 FITTING REDDENING LAWS 

We will use a similar method to A14 to fit different 
parametrised reddening laws to the measured colour ex- 


E{X -V) = {X-V)-{X- V)^ , 


( 2 ) 


where {X — 1/)q is the assumed X — V colour for the unred¬ 
dened objects. 

In this work we will study different extinction laws and 
how well they describe the observed SN la reddening. Each 
extinction law, A(A;a), depends on a set of parameters a 
(e.g. E{B— V) and Rv), that will be fitted to the observed 
colour excesses, E{X — V)p, by minimising 


2 _ [E'(A V)p (Axp ^Vp)] 

X n Vp 


(3) 


Here Axp and Avp are the predicted extinctions in the X 
and V filters for phase p and can be calculated as 


Axp = -2.5 log 1 


/ Tx{X) ■ . SoiX-,p)XdX 

jTx{X)So{X-,p)XdX 


(4) 

if the effective filter transmission, Tx{X), and the spec¬ 
tral energy distribution (SED), 5o(A;p), of the unreddened 
source are assumed to be known. Further, (TXpVp s-re the 
uncertainties added in quadrature. This includes the mea¬ 
surement errors, but will in most cases be dominated by the 
intrinsic colour uncertainties discussed in § |4.1| 

In addition to the extinction law parameters a, we will 
also minimise Equation (|^ with respect to the parameters 
and SB that together determine the phase, p. The ob¬ 
served values of these parameters are extinction dependent 
(see e.g. Leibundgntj 1988 Phillips 1993 Nugent, Kim fc 
Perlmutter|2002| A14) since broadband measurement of an 
object suffering from extinction will effectively probe redder 
wavelengths than the observation of the same unreddened 
object would have. The SN la light curve decline rate varies 
with wavelength, and slower decline rates are in general ex¬ 
pected to be observed for reddened objects. 

Similarly, the Galactic extinction, A^"^, as well as the 
Kx corrections, are also properties measured through broad¬ 
band filters and therefore also depend on the SED of the 
source, as seen from e.g. Equation Since the observed 
SED in the Milky Way is affected by the extinction in the 
SN host galaxy, these properties will also depend on the 
reddening law. The effect will be particularly significant in 
the UV, where the wavelength dependence for the extinc¬ 
tion laws is steep. For example, the Galactic extinction in 
the Swift/ uvm2 band for the line of sight to SN 2014J would 
have been AA[^ « 0.3 for an unreddened source while it is 
A^'^ « 0.1 for SN2014J. 

We take this effect into account by fitting each redden¬ 
ing law, A(A; a), iteratively and update A^'^ and Kx using 
the fit a values in each iteration. That is, after each itera¬ 
tion, we re-calculate the observed colour and colour excesses 
using Equations 0 and 0 before re-fitting the reddening 
law using Equation ([^ . The procedure is repeated until the 


change in the fitted parameters is less than 1 % between 
iterations. 

As mentioned above we use the average value of Rv ~ 
3.1 when correcting for the Galactic extinction. For the 
majority of the SNe the Galactic extinction is negligible 
compared to the host reddening, with the exception of 
SN 2012et. For this SN we also tried carrying out the fits 
using the extreme values of Rv = 2.2 and Rv ~ 5.8 ob¬ 
served in the Milky Way (e.g. COM, Fitzpatrick||1999 1 and 
concluded that this could impact the fit values with up to 
one statistical standard deviation. 


4.1 SN la intrinsic colours and SED 

In order to minimise Equation ^ , we need to make assump¬ 
tions of the intrinsic colours, {X — V)^, and SED, Sq(X\p) 
of the unreddened objects. The reddening laws for SNe la 
are typically derived by either comparing them to individual 
objects that show similarities in light curve properties and 
spectral evolution (see e.g. Krisciunas et al.]|2006[ A14), or 
to SN la colour and SED templates. 

In Figure[^the H07 template is plotted (dashed red) for 
the optical colours. Although the H07 template does extend 
from the UV to NIR, the data it is based on are sparse at 
the endpoints of this range. N08 studied the intrinsic optical 
colours between —10 and -|-50 days from R-band maximum. 
Two of their colour laws are shown (in blue) for normal, 
Sb = 1, SNe la in Figure together with the colour dis¬ 
persions they derive. Although these templates provide ex¬ 
cellent coverage at optical wavelengths none of them cover, 
to high accuracy, the full wavelength range required for the 
UV-NIR analysis in this work. 

The best studied unreddened SN la to date is SN 2011fe, 
discovered by the Palomar Transient Factory in the nearby 
spiral galaxy MlOl (PTFllkly, Nugent et al.|20lT l. Its close 
proximity allowed detailed spectroscopic and photometric 
observations over a broad wavelength range from the UV 
(Brown et al.||2012 Mazzali et al.|[2C)14l, through the opti¬ 


cal (e.g. Munari et al.||20T3 l to the near- (Matheson et al. 
|2012[ [Hsiao et al.|2013| , and mid- ( [McClelland et al.|2013p,' 


infrared. The SN was also targeted in the far-IR (Johans- 
son, Amanullah, fc Goobar|2013 1 and radio (Chomiuk et al. 


2012b|, but was not detected at these wavelengths. 


SN2011fe is a normal SN la in both the optical and 


the NIR (Munari et al. 2013 Matheson et al. 20121 and 
does not show any spectroscopic peculiarities ([Pereira et al.[ 


20131. Further, the low Galactic and host galaxy reddening 


along the line-of-sight, E{B — 1/)mw = 0.011 ± 0.002 mag 
and E(B — V)host = 0.014 ± 0.002 mag (deduced from the 
integrated equivalent widths of the Na l D lines; [Patat et alT] 
2013[|, makes it an excellent comparison object for studying 


reddening of SNe la. The colours of SN 20life shown in Fig- 
ure[^(solid black lines) have been obtained by combining the 
available UV-NIR data, as described in Appendix [B[ based 
on the measurements (black triangles). We have also com¬ 
bined the spectroscopic data to create a daily sampled SED, 
shown in Figure [Bl| From the excellent B — V match be¬ 
tween SN 2011fe and H07 we can also conclude that any po¬ 
tential discrepancies introduced by using the SNooPy model 
for the U-band is not likely to have any significant effect on 
the derived reddening laws when the colours are compared 
to the corresponding SN 201 Ife colours. 


© 2015 RAS, MNRAS 000, 



































































Diversity in extinction laws for SNe la 11 


In order to properly compare the colours of reddened 
SNe la to SN 2011fe, and fit reddening laws, we also need 
accurate estimates of the expected intrinsic colour disper¬ 
sions. It is desirable to take into account how these vary in 
time, and how they are correlated both in time and with dif¬ 
ferent colours. However, since we lack this information for 
all colours involved in the analysis, we adopt a simplified 
approach where all colours are treated equally with phase- 
independent dispersions for each colour. We further follow 
the procedure from A14 and assume that the colour uncer¬ 
tainties between SN phases are completely correlated and 
that there is no correlation between different colours. 


Milne et al. (|2013 M13 from hereon) following the work 


of Brown et al.| (20101 and Milne et al. (20101, studied 
the UV—optical colours of 23 SNe la observed by Swift. 
They reported a dispersion of ~ 0.3 mag in UV—n for their 
E{B — V) < 0.2 mag sample after correcting the data for 
extinction and the result does not seem to be affected sig¬ 
nificantly by the choice oi Rv- We adopt this dispersion, 
shown as grey bands around the SN 2011fe colour template 
in Figurej^ for the uvm2 — V, F225W— V, and F275W— V 
colours. For the F336W — V and U — V we adopt a disper¬ 
sion of 0.1 mag based on the U — V dispersion from NOS. 
Note that since these dispersions were derived for unred¬ 
dened objects, and the dispersion of SNe la colours decrease 
with wavelength, they can be considered as conservative for 
highly reddened SNe. When reddened SNe la are studied 
using the WFC3/UVIS filters, the observations will effec¬ 
tively probe redder wavelengths compared to similar studies 
of unreddened objects as illustrated in e.g. Figure [dTJ and 
as a consequence the intrinsic dispersion of the former obser¬ 
vations can be expected to be lower than the latter. On the 
other hand, even larger dispersions have been also observed 
in the UV, and we will use these intrinsic uncertainties un¬ 
der the assumption that all the SNe we study are similar 
to the reference SN 2011fe, which we will discuss further in 
§ 1 ^ 

B14 derive a dispersion for the pseudo-colour Bmax — 
Umax of 0.06 mag using a Cauchy Prior for the colour distri¬ 
bution of their SNe la observed with the Carnegie Supernova 
Project. This is consistent with the phase-dependent results 
of NOS who derive a. B — V dispersion of < 0.1 mag, which 
we use here, within the range —10 to -1-35 days. For the 

V — R and V — i colours we adopt the NOS dispersions of 
O.OS mag. 

The NIR colours of SNe la typically show low dispersion 
(see e.g. Meikle 20001 and for the NIR-optical V — J and 

V — H colours we can estimate the dispersion in our phase 
range by using all SNe with E{B — V) < 0.2 mag from the 
’’Gold Sample” in Stanishev et al. (20151 and find 0.12 mag 
and 0.13 mag respectively. For V — Ks we conservatively 
adopt 0.2 mag based on the studies carried out at maximum 


(B14) and the dispersion of the /fs-band (Friedman 2012 
[Friedman et al.]|2014[ ). 

The observed colour excesses, E{X — V), for each red¬ 
dened SN, together with the adopted colour dispersions, 
shown as background regions, are plotted in Figure 
SN 2012bl has not been plotted in the figure. It will be omit¬ 
ted from the extinction analysis due to its lack of reddening. 

If the colour excess is observed at single wavelengths, 
i.e. E{Xi — A 2 ) = — A\ 2 , it will be time-independent for 

interstellar dust extinction. However, for broad-band pho¬ 


tometry, we can expect E{X — V) to vary with time which is 
illustrated by e.g. the UV colours in Figure]^ As the intrin¬ 
sic colour of the source varies with time, the effective wave¬ 
lengths of all broad-band filters will change as well. Further, 
since the first term on the right-hand side in Equation § is 
also affected by a reddening law, any change in the effective 
wavelength of the filters X and V will affect this term more 
than the second term, which will induce a time-variability for 
E{X — V). For most filters, the time-variability will be neg¬ 
ligible except in the UV (where the extinction has a steep 
wavelength dependence), and in particular for filters with 
red tails, e.g. F225W (see Appendix]^, the effect is signif¬ 
icant. This is illustrated by the bottom panel of Figure [DT] 
where the effective wavelength of three WFC3/UVIS UV fil¬ 
ters have been plotted for a SN la at maximum for different 
extinctions. 


4.2 Reddening laws 


We test three different extinction law parametrisations. In 
addition to the widely used CCM law modified by jQ’Donnellj 
( |1994[ hereafter CCM-l-0), which has been derived from 
studying different lines-of-sight in the Milky Way, we also 
fit the parametrisation from F99 that has been derived in 
a similar manner. For each of these we fit both the colour 
excess E{B— V), which relates to the optical depth, and the 
ratio of the total-to-selective extinction, Rv- In the Milky 
Way the value of Rv typically varies between Rv ~ 2.2- 
5.8, for different lines-of-sight (e.g. CCM, Fitzpatrick|1999 1, 
despite this we allow extrapolations of Rv within the range 
Rv = 0.5-8 when we fit the laws to the SN colours. 

We also test a simple power-law model ( Goobar|2008 1, 


A(A; Av, a, f)) = Av 


1 — a -1- a I -— 
Av 


where the reference wavelength. Ay, will be chosen as Ay = 
0.55 /rm. The CCM-I-O and F99 laws can be approximated 
by a power-law in the optical and NIR range, and jCoobarj 
(20081 showed that the observed reddening law of an object 


embedded in circumstellar dust can be approximated by the 
expression above with a — 0.8 (a = 0.9) and P = —2.5 
(/3 = —1.5) for Milky Way (LMC) like dust. The reason 
why the original extinction laws are not preserved in the CS 
scenario is due to the geometry where multiply scattered 
photons on the CS dust will reach the observer while this 
is extremely unlikely for an ordinary interstellar dust geom¬ 


etry. Note that the values in Goobar (20081 were obtained 


under the assumption that the light source is constant in 
flux and colour. For time-dependent light sources the ob¬ 
served reddening is expected to vary in time which has been 
investigated for SNe la by e.g. jAmanullah fc Goobar'| ( |2011| ). 

When the power-law is fitted to colour excesses, E{X — 
V) = Ax — Av, Av and a are almost completely degenerate, 
and we are in fact only sensitive to the product a ■ Av- We 
will for this reason fix a to a = 1 in all fits below, which also 
has the implication that the resulting Av can no longer be 
interpreted as solemnly being the extinction in the U-band. 
We have further checked that fixing a does not affect the 
fitted values of fi, and that we obtain consistent values of 
the product a ■ Av when a is allowed as a free parameter, 
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Figure 6. The observed colour excess, E(X — V), against phase for all SNe and filters. The predicted colour excess from the SN2011fe 
template and the fitted F99 law, presented in Table is shown with solid lines and the regions surrounding them mark the adopted 
intrinsic colour uncertainties. The uncertainties are also indicated by open brackets for the three SNe with the least reddening. The 
phase, p = — was calculating based on the fitted and sb values from Tablej^ The apparent E{X — V) time-dependence 

of the F225W filter for e.g. SN2014J originates from the red tails of these filters which causes the observed colour excess to vary with 
the intrinsic colour of the source. 
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although the individual values of a and Av are in this case 
of course poorly constrained. 

In A14, the SALT2 ( Guy et al.|2007 1 law from Betoule 


et al. (20141 was also studied. This is a colour law and not 


an extincton law, and is further only defined for wavelengths 
A < 8000 A. It can therefore not be compared on an equal 
footing to the other laws and will not be studied in detail in 
this work. We will however compare it to the reddening or 
our low extinction SNe in §|6.5| 


4.3 Results 

The fitted parameters for the three extinction laws are pre¬ 
sented in Table [5] for each individual SN. The extinction 
laws were fitted to all available colours between phases —10 
and -1-35 days and the predicted colour excesses using the 
best fitted F99 laws are also shown in Figurej^ The revised 
photometry of SN 2014J and updated intrinsic colour model 
with respect to A14 only had a minor impact on the fitted 
results for this highly reddened SN, and the values presented 
here are consistent with A14. 

Further, all three reddening laws capture the general 
wavelength dependence of the observed colour excesses with 
the possible exception of SN 2012cg which will be discussed 
below. This shows that the observed SN colours from UV 
to NIR can indeed be described by the colours of SN 201 Ife, 
with the adopted intrinsic dispersions, together with an ex¬ 
tinction law with only two free parameters. 

The six SNe span a broad range both in reddening, 
E{B — V) = 0.2-1.4 mag, and in Rv ~ 1.4-3.0 which con¬ 
firms the findings from previous optical and NIR studies 
(e.g. B14), that point to a diversity of observed reddening 
laws of SNe la. The significant difference in the derived ex¬ 
tinction laws is also illustrated in Figurej^where the relative 
extinction, A\/Av, for the F99 law has been plotted for two 
different values of Rv- The diversity is particularly striking 
blue-wards of the [/-band (A“^ > 2.5 /im“^), emphasising 
the power of using UV data to study diversity in extinction. 

In the Figure we also show the relative extinction, 
Ax/Ay = 1-1- E{X — V)/Av, for the four SNe with 
E{B — V) > 0.2 mag against inverse wavelength. Here, 
E{X—V)\s calculated as the weighted average for all phases 
and Av was obtained from the fitted F99 parameters as 
Av ~ Rv ■ E{B — V). Because of this, the derived values 
for Ax!Ay are model-dependent and will be different for 
the different extinction laws. Further, since the extinction 
Ax is obtained from comparing observations in the filter X 
of reddened and unreddened sources, and the measurements 
correspond to different effective wavelengths, it is non-trivial 
to derive a specific wavelength corresponding to Ax/Ay - For 
these reasons, we follow the procedure from A14 and plot the 
AxjAv points in a manner where we chose to preserve the 
residuals between the fitted model and the data points from 
Figure]^ This is obtained by plotting each point, AxjAy, 
at the wavelength where the residual to the extinction law 
is the same as the corresponding weighted average residual 
of E[X — V) over all epochs. We emphasise that this plot 
is only a way of visualising the data and the reddening laws 
and that all extinction parameters were obtained by min¬ 
imising Equation i.e. not by fitting the extinction laws 
to the data points in Figure 

One observation that can be made by comparing the re- 



Wavelength ^ (^m 

Figure 8. Similar to Figure but here the power-law fits from 
Table[^are shown with the observed reddening. The wavelengths 
of the data points will differ from Figure]^ due to the way these 
plots are created, as explained in the text. Note that all fitted 
parameters described in this work were obtained by minimising 
Equation §, and not by fitting the extinction laws to the data 
points in this plot. 


suits in Table and Figure is that although all extinction 
laws appear to be general enough to describe all SNe, the 
F99 law shows tension for SNe 2012bm and 2014J. These 
SNe have in common that they both show a fair amount of 
reddening {E{B — V) > 0.5 mag). The more dust along the 
line-of-sight to a SN, the more sensitive the fit will be to 
the extinction law, and less to possible intrinsic SN colour 
variations. In other words, we can expect any deviation be¬ 
tween the assumed parametrisation and the actual observed 
extinction law to become more significant the redder the 
SNe are. It has already been pointed out that since the F99 
law is empirical in nature, it may have features that do not 
properly extrapolate to low values of i?i/. This could be the 
case for SN 2014J, where the F99 law with Rv = 1.4 leads 
to a more prominent feature around ~ 1.5 /im“^ (~ 6700 A) 
than what is supported by the data. In Figure|^we compare 
the power-law fits for the four most reddened SNe, and here 
we can see how this is consistent with all the data points 
over the full wavelength range for SN 2014J. 

From the reduced in Table we can conclude that 
the power-law model provides better fits than F99 for both 
low-i?v SNe, 2012et and 2014J, while the opposite is true for 
the Rv < 3 SNe 2012cg, 2012cp and 2012cu. In particular, 
the F99 and CCM+O law seem to provide nearly perfect 
fits to the observations of SN 2012cu, which is also seen in 
FigurejT] SN2012bm shows tension for all fitted laws. 

The CCM-fO law performs better in terms of goodness- 
of-fit for several cases, e.g. for SN 2014J. However, this law 
does not extrapolate well to low Rv values which give rise 
to a sharp, discontinuous, knee when the NIR and optical 
parametrisations of the law are merged. 
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2012cg 

2012et 

2012cp 

2012bm 

2012CU 

2014J 


E{B- V) 

0.15(0.02) 

0.17(0.03) 

0.35 (0.03) 

0.46(0.04) 

0.99 (0.03) 

1.36(0.02) 

Oi 

cni 

+ 

Ry 

Q iy+0.9 
^•^-0.7 

1 rr-|-0.6 

-*-•'-0.5 

3 

'^•0'-0.4 

3 q+0-5 

d d 
+ 1 
00 

-*-•^-0.1 

^max 

56080.6(0.3) 

56191.0(0.5) 

56081.1 (0.3) 

56018.0(1.9) 

56104.5(0.4) 

56687.9(0.2) 

Sr 

1.12(0.02) 

1.01 (0.03) 

1.18(0.02) 

1.24(0.12) 

1.04 (0.03) 

1.14(0.01) 



1.67 

1.50 

1.60 

2.83 

0.63 

2.38 


E(B- V) 

0.15(0.02) 

0.17(0.03) 

0.36 (0.03) 

0.48(0.04) 

1.00 (0.03) 

1.40(0.02) 

o 

+ 

S 

Ry 

^•'^-0.7 

1 4+0-8 
-*-•^-0.7 

CO 

b 

1 + 

o o 

2-8+.1 

rs q + 0.2 

1 

^•^-0.1 

xmax 

56080.6(0.3) 

56191.1 (0.5) 

56081.1 (0.3) 

56019.1 (1.7) 

56104.7(0.4) 

56688.1 (0.2) 

o 

o 

SB 

1.12(0.02) 

1.01 (0.03) 

1.18(0.02) 

1.18(0.10) 

1.03 (0.03) 

1.13(0.01) 



1.60 

1.34 

1.42 

2.14 

0.91 

1.42 


Av 

0.44(0.19) 

0.30(0.13) 

1.51 (0.40) 

1.47(0.43) 

3.59 (0.29) 

1.97(0.10) 

1 

/3 

-1.2 (0.4) 

-1.7 (0.5) 

-0.8 (0.2) 

-1.2 (0.4) 

-1.0 (0.1) 

-2.0(0.!) 


amax 

^B 

56080.6(0.3) 

56191.1 (0.5) 

56081.2 (0.3) 

56020.4(2.7) 

56104.8(0.4) 

56688.1 (0.2) 

0 

SB 

1.12(0.02) 

1.00(0.03) 

1.18(0.02) 

1.16(0.15) 

1.02(0.03) 

1.13(0.01) 


yjv 

1.79 

1.17 

1.95 

2.13 

1.78 

1.65 


Table 5. The best fitted parameters to all measured colours between phases —10 and +35 days for the reddening laws investigated in 
this work. For the power-law we have fixed a to a = 1 to break the degeneracy between this parameter and Ay. For each value we 
quote the 68 % level fitted uncertainty when the parameters are considered individually. The degrees of freedoms were calculated under 
assumptions that the colours are independent measurements, which is not the case. 


Ks H J i R V B U F275W 

F336W F225W 



Wavelength ^ {^m M 

Figure 7. The normalised extinction, Ax/Ay for the four SNe with E{B — V) > 0.2 mag. Here, Ax/Ay = 1 + E{X — V)/Ay where 
E(X — V) is calculated as the weighted (with the measurement uncertainties) averages for all epochs, and Ay is obtained from the 
fitted F99 parameters in Table[^ Following A14, each data point has been plotted at the wavelength where the residual to the extinction 
law is matching the corresponding residual between E{X — V) and the fitted reddening law in Figure]^ Note that all fitted parameters 
described in this work were obtained by minimising Equation §, and not by fitting the extinction laws to the data points in this plot. 
The filter passbands for the filters used for the minimisation are shown at the top. 
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SN 

Phase 

E{B- V) 

Rv 


SN2012cg 

[-5,+5] 

0.16 

(0.02) 

1 y+O-S 
-*^•'-0.6 

1.51 

[+10,+20] 

0.13 

(0.03) 

'^•'^-1.0 

1.34 

SN2012cp 

[-5,+5] 

0.33 

(0.03) 

1 q+0.6 

1*^-0.5 

1.63 

[+10,+20] 

0.34 

(0.03) 

9 9+0.8 
^•^-0.6 

0.55 

SN2014J 

[-5,+5] 

1.35 

(0.03) 

1 

-*^•^-0.1 

1.81 


[+10,+20] 

1.32 

(0.03) 

-‘^•'^-0.1 

1.56 


Table 6. Fitted F99 parameters using only data between —5- 
+5 and +10—1-20 days from maximum, respectively. While SNe 
2012cg and 2014J both have UV-NIR observations for both in¬ 
tervals, the fits to SN 2012cp are only based on UV-optical data. 
The parameters and sg were fixed to the values given in 

Table 


5 SEARCHING FOR CIRCUMSTELLAR DUST 


All the SNe studied here show a colour excess evolution 
that is within the adopted dispersion with the exception 
of SN 2012cg. This SN reveals a minimum in its NIR colour 
evolution around two weeks past max which is possibly ac¬ 
counted to the diversity observed in NIR SN la light curves 
between the first and second peaks. This is within the as¬ 
sumed intrinsic uncertainty although the assumption that 
these are fully correlated between epochs does not hold in 
this case. What is perhaps even more striking are the ~ 2a 
variations of the uvm2— V evolution which suggests that this 
colour evolve slower than the corresponding 2011fe colour, 
resulting in an apparent decrease with time of the colour 
excesses. 

The fact that the observed colour evolution of 
SN 2012cg is not consistent with a single extinction law is 
also illustrated by Table where the F99 law has been fit¬ 
ted for two phase intervals: around maximum and between 
10 and 20 days past The fits were carried out for the 
three SNe that had UV data in both these intervals while 
keeping the fitted values of and sg fixed to the values 
given in Table The fitted parameters of SNe 2012cp and 
2014J are within errors between the two epochs while this 
is not the case for SN 2012cg. 

Evolving colour excesses, or a reddening law that 
changes with time, is predicted if circumstellar dust is 


present in the SN environment. Amanullah & Goobar (20111 


studied this effect at optical wavelengths and we can ex¬ 
tend that analysis, using the same tools, to also cover UV 
and NIR in order to test the observed colour evolution of 
SN 2012cg against the expectations for a CS scenario. 

In Figure the same colour excesses of SN 2012cg that 
are plotted in the upper left panel of Figure are shown. 
The solid lines represent the result when we fit a combined 
F99 and CS dust law to the bluest colours, uvm2 — V and 
F225W — V, while keeping and sg fixed to the values 
from Tablej^ The best fit values for the F99 law are E{B — 
V) ~ 0.2 mag and Rv ~ 3.0. The CS dust in this case 
consists of dust modeled to match the average extinction 
properties in the Milky Way ( Draine||2003 1 and distributed 



Figure 9. The measured colour excesses for different phases for 
SN2012cg are taken from Figure]^ Overplotted are the best fits 
of a combined circumstellar and F99 dust law, where the former 
has been adopted from [Amanullah &: Goobar] ) |2011| l where the 
dust is distributed in a thin spherical shell. The solid lines shows 
the best fit when only the uvm2 — V and F225W — V colours are 
considered, while the dashed lines show the fit to all colours. The 
former is obtained with a CS dust shell at radius r^ust ~ 10^^ cm 
from the SN that give rise to an extinction of Ay ~ 0.2, while 
latter fit is consistent with no CS dust. 


in a thin spherical shell. Here the best fit yields an extinction 
of « 0.2 for a shell at a radius of ruust ~ 10^^ cm. 

Although adding Milky Way-like CS dust to the picture 
can give an explanation for the time-evolution of the bluest 
colours, it does not provide an improvement beyond a sin¬ 
gle F99 law when all colours are considered simultaneously, 
which is shown by the dashed lines in Figurej^ Here the best 
fit is obtained for « 0.0. The explanation for this is 

that the specific CS dust model we test also predicts evolu¬ 
tion of the redder colours which in this case is not supported 
by the data. If CS dust is present, the total extinction. Ax, 
can only decrease with time due to late-arriving scattered 
photons. From the results in Table[^ we see that this is not 
the case for Ay, although the value of Ay is in this case 
model dependent. 

It is also possible to detect CS dust by studying the 
time-evolution of multi-epoch high-resolution spectroscopy. 
This was mentioned in §[^and is the observational signature 
that has been used to claim detections of CS dust around 
SNe la. 

For SN 2012cg we obtained high-resolution FIES spec¬ 
tra for three epochs which are listed in Table and corre¬ 
spond to phases —7, 1 and 13 days with respect to All 
spectra contain well resolved unsaturated Na l D features 
(Figure]^, which are shown in Figure 10 The bulk of the 
distinguishable features are likely part of the visible inter¬ 
stellar medium (ISM), but it can not be excluded that Na l 
associated with CS medium could contribute to the observed 


© 2015 RAS, MNRAS 000,[iH33| 

























16 


R. Amanullah et al. 


Velocity (km/s) 

350 400 450 500 550 600 



Figure 10. Na I D absorption doublet of SN 2012cg. The rest frame wavelength of Na I D in NGC 4424 (solid lines) and for the 
approximate stellar velocity along the line-of-sight (dashed) are indicated. The dotted vertical line indicates a sky emission line. In 
Section [^the profile of Na I D is discussed in detail. 


profile. A full description of the Na l D profile and summary 
of other interstellar absorption features is presented § |6.4[ 

The FIES spectra were normalised by fitting 3rd-order 
polynomials to the continua bracketing Na l D. No telluric 
corrections were performed on the spectra, due to the lack 
of standard star observations on some of the epochs. We fit¬ 
ted Gaussian profiles to the weighted average of the three 
epochs, and find that the Na l D line ratios (D1/D2) range 
from 2-1.2 for individual features. The deepest features have 
low ratios, indicating that they are not optically thin. Based 
on a standard star spectrum taken on the second epoch, we 
identify telluric features overlapping with Na l D2, where 
as Na I D1 appears to be located at a less affected part of 
the spectrum. We therefore focus most of the time-variation 
analysis on Na l Dl, because it is less contaminated by tel¬ 
luric features and it is more likely to be optically thin due 
to the lower absorption cross-section of this transition. We 
measure the total equivalent width of Na l Dl of the respec¬ 
tive epochs to be, 702 ± 9, 711 ± 32 and 685 ± 20 mA. 


Since no significant time-variations of Na l D absorp¬ 
tion between different epochs could be detected beyond the 
noise level, we attempt to use the non-detections to con¬ 
strain possible CS dust models assuming a thin sphere of 
CS material shell distribution around the SN. Using the Na l 
photo-ionisation model described by [Borkowski, Blondin fc 


Reynold^ 1 2009[) and recently applied to SN 2014J in |Gra- 


ham et al. (2015bI, we can exclude CS material, and thus 


also CS dust, at certain radii from the SN. If the Na l is 
optically thin, a change in equivalent width of Na l D is di¬ 
rectly proportional to a decrease in column density. Further¬ 
more, the product of the integrated photon count flux above 
the ionisation energy of Na l (2400 A) and its correspond¬ 
ing ionisation cross-section is proportional to the fractional 
decrease in column density. We can thus compute the frac¬ 
tional decrease in equivalent width of Na l D expected due to 
photo-ionisation of CS material at a given radius from a SN. 
In this model, CS Na l can be excluded which is far enough 
away from the SN to not have been fully ionised by the 
first epoch and close enough to see changes by the following 
epochs. The diversity of SNe la spectra in UV implies that 
the results will be sensitive to the spectral templates used. 
Noting that the STIS spectra of SN2012cg resemble those 


of SN2011fe, we use the SN2011fe spectral template de¬ 
scribed in Appendix]^ However, the H07 template extends 
to shorter wavelengths and earlier times than the SN 201 Ife 
template. We therefore extrapolate our spectral template by 
H07, whereby the flux is rescaled for continuity. This slightly 
increases the photo-ionising flux. 

In Figure El the model fractional decrease of Na l at 
different radii from a SN are shown. We determined the 
inner radius of excluded CS Na l of 2012cg by assuming 
that an absorption feature deeper than three times the root- 
mean-square of the noise would be detectable in the ear¬ 
liest spectrum. This yields a distance of rdust < 2 • 10^® 
cm within which the Na l column density would have de¬ 
creased below detection levels since explosion. With respect 
to the first epoch, the fractional change of the equivalent 
width of Na l Dl is considered to set an upper limit on the 
change in column density. Within du errors, the last epoch 
excludes the models in which the dust is located closer than 
J’dust > 3T0^® cm from SN 2012cg. Using the total equivalent 
width of Na l Dl of course has the caveat that it considers 
the entire measured column density to be situated at one 
radius from the SN, which we know cannot be true based 
on the profile. We therefore also consider how much a sin¬ 
gle feature must change to be detected given the signal-to- 
noise of the spectra. The limits obtained by the absence of 
change beyond three times the root-mean-square noise are 
comparable to those set by the fractional change of the total 
equivalent width. Lastly we note that the limits are not valid 
if there is CS material that is not optically thin, since the 
equivalent width then does not change linearly with column 
density. 


In summary, although the existence of Milky Way-like 
CS dust at rdust ~ 10^^ cm is consistent with the ob¬ 
served time-evolution of E{uvm2 — V) and non-detection of 
time-varying Na l D absorption for SN 2012cg, the minimal 
E{X — V) evolution of the redder wavelengths disfavours 
this explanation. This suggests that CS dust can only be 


present either in very small amounts, Ay 


< 0.05 or 


if it is located at larger radii rdust ~ 10^® cm for which 
the time-delays are much larger than the time-scales stud¬ 
ied here and are not expected to affect the observed colours 
(Amanullah fc Goobar|2011 1. Further, non-detections in ra- 
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Phase (days) 


Figure 11. Na i ionisation models for thin shells of CS mate¬ 
rial at different distances from a SN, suggesting the absence of 
circumstellar material at distances between r^ust = 2 • 10^^ - 
3 • 10^^ cm around SN2012cg. The upper panel shows the frac¬ 
tional decrease of Na i due to ionization as a function of time at 
radii ranging from rdust = 10^^ — 3-10^® cm. The arrow indicates 
an approximate upper limit for Na I around 2012cg to be ionised 
beyond detection level by —8 days from B maximum based on 
the S/N around the Na i D doublet of the first FIES spectrum. 
The lower panel shows the fractional change of Na i with respect 
to the first epoch. The fractional change of equivalent width of 
Na I D1 with 3 (j error bars set an approximate lower radius limit 
at which CS material could be present while the ionisation effects 
remain undetected. 


dio (Chomiuk et aL|2012al and far-IR (Johansson, Amanul- 


|lah fc Goobar|2013^ are other observational signatures sug¬ 
gesting that SN 2012cg exploded in an environment free from 
CS material. Here the far-IR observations are sensitive to 
pre-existing CS dust that would be heated by the explosion. 


6 DISCUSSION 

6.1 Using SN 2011fe as a colour template 

One possible explanation for the time-evolution of the 
SN 2012cg UV — V colour excess could be that it is intrinsi¬ 
cally different from SN 201 Ife. We have already argued that 
all the SNe in our sample, are normal SNe la by comparing 
their overall spectra (see Figure]^ and their SNID classifi¬ 
cation, but we can also study individual features which has 
proven to be a useful approach for sub-classifying normal 
SNe la. 

The properties of individual features can be quanti- 


SN 

Phase 

(days) 

pWl 

(A) 

pW6 

(A) 

pW7 

(A) 

Type 

2011fe 

+0 

94(1) 

14(1) 

94(1) 

GN 

2012bl 

-1 

110(2) 

1(1) 

89(2) 

CN 

2012bm 

+9 

125(2) 

13(1) 

76(2) 

CN 

2012cg 

+2 

103 (4) 

9(3) 

75 (4) 

CN 

2012cp 

+2 

98(1) 

10(1) 

74(1) 

CN 

2012CU 

+9 

70(1) 

24(1) 

99(1) 

CN 

2012et 

+2 

123(2) 

11(1) 

133(2) 

BL 

2014J 

-1 

155(2) 

14(1) 

100 (3) 

CN 


Table 7. Calculated pseudo-equivalent widths {pW) based on 
the optical spectra. The p W naming convention from |GaravinI] 
|et al.| l (2007[ l is followed with pWl corresponding to Ca II H&K, 
and pWd and pWJ corresponding to Si II A5972 A and Si II 
A6355 A respectively. The quoted uncertainties are the statis¬ 
tical errors. Systematic uncertainties from estimating the pseudo 
continuum could be dominating in some cases. The classification 
as ’’core normal” (CN) or ”broad-line” (BL) is based on the cri¬ 
teria from [Blondin et al.| l |2012[ | . Given that these criteria were 
derived based on spectra around maximum, the classification of 
SN 2012cu should be considered less certain. The measurement 
of SN 2014J was carried out on a spectrum published in|Goobar| 
|et al.| | |20l4 . See the text for details. 


fled by e.g. measuring their velocities and pseudo-equivalent 
widths (pW), where pseudo refers to the fact that the equiv¬ 
alent widths are obtained using a pseudo continuum since 
proper continua are absent in SN la spectra. The pseudo 
continuum can be defined with a straight line between the 
two flux peaks surrounding an absorption feature and the 
pW IS then calculated as the integral of the spectrum flux 
relative to the continuum. Using this method, the error of 
the measurement will typically be dominated by the sys¬ 
tematic uncertainty introduced by the determination of the 
pseudo-continuum. 

[Branch et al.| ( |2006| ) suggested a classification scheme 
based on pseudo-equivalent width measurements of the ab¬ 
sorption features near 5750 A and 6100 A (following the 
naming convention from [Garavini et al. j pOOTj ) we refer to 
these as pWG and pWl respectively), that can be associ¬ 
ated with Si II A5972 A and Si ii A6355 A as shown in 
Figure [^ They identify four different groups when consider¬ 
ing pWQ and pWT for their sample where in particular the 
’’core normal” (CN) SNe, to which SN2011fe belongs, are 
tightly clustered and show a high degree of general spectral 
homogeneity. 

In Table [T] we show measurements of p W6 and p W7 
for spectra close to maximum after correcting our spectra 
for host galaxy reddening (although this is not expected to 
have any major impact on the results as discussed injNordinj 


et al. 20111 using the fitted reddening laws from §4.3 


We 


also presen t the p W for Ca ll H&K (p Wl following Garavini 


et al. 


20071 and note that all SNe have pWl < 150 A. It has 


been shown that pWl correlates with the intrinsic U — B 
colour (Folatelli et al. 2013 Burns et al. 20141 where SNe 
with pWl> 150 A are intrinsically redder than the bulk of 
normal SN la. 

Using the classification criteria from [Blondin et al.j 
(20121 we can classify all SNe as CN except for SN 2012et 
which lands in range that Branch et al. {20061 defines as the 
”broad-line” (BL) group. These SNe la have broader and 
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deeper 6100 A absorption but are in most aspects not very 
different from CN SNe. When the observed spectra from the 
two groups are compared with synthetic spectra generated 
with SYNOW ( Branch et al.|2003 l similar photospheric ve¬ 
locities and excitation temperatures can be used to describe 
both the CN and BL groups ( Branch et al.|2006 |. 

Given these results, there is no evidence for any signifi¬ 
cant discrepancy between the optical properties of our SNe 
and SN 2011fe, with the exception for the additional source 
of luminosity observed in the very early optical lightcurve 
of SN 2014J (Goobar et al. 20151. However, several stud¬ 


ies have shown that normal SNe in the optical still show 
dispersion in the UV and M13 argued that their sample 
could be divided in up to four groups. Here, SN2011fe be¬ 
longs to what they define as the ”near-UV (NUV) blue” 
group which they conclude is on average 0.44 mag bluer in 
u — V than the reddest, ”NUV-red”, group. They consider 
SNe with E{B — V) < 0.25 mag, and this constraint could 
still allow extinction to contribute with up to 0.9-1.3 mag 
(depending on the extinction law) to the uvm2 — v colour, 
which may offer an explanation for the significant colour dis¬ 
persion. However, M13 show that the two groups are distin¬ 
guishable even after the individual SNe have been corrected 
for reddening. 

Further, they also conclude that the groups, on av¬ 
erage, show different spectroscopic properties. All NUV- 
blue SNe have ’’normal” photospheric expansion velocities 
around maximum when quantified based on the Si ll A6355 A 
feature, while the remaining groups are populated with both 
’’normal” (NV) and ’’high” (HV) veloc ity SNe. They use 
the definition from Wang et al. (20091, where SNe with 


Si II velocities of usi n > 11,800 kms“ at maximum are 
classified as HV and the remaining as NV. Although there 
may not be a sharp distinction between the two classes (see 


e.g. Silverman, Kong & Filippenko 20121, the M13 results 


are consistent with previous studies showing that an opti¬ 
cal colour-velocity dependence exists with HV SNe being on 


average intrinsically redder ( 

Foley & Kasen 

20111 |Blondin 

et al.[[2012| [Folatelli et al.[[2013[ [Mandel, Fo 

ey & Kirshner 


2014). On the other hand, exceptions from these findings is, 
SN 2011de, the brightest SN la observed ( Brown|2014 1, and 
the results by Wang et al. (2012), showing, using the HST 


Advanced Camera for Surveys, that the UV-optical colours 
of the HV object SN 2004dt are significantly bluer than the 
NV SN 2005cf. 

We have measured the Si ll velocities, which are pre¬ 
sented in Figure for all spectra we have obtained of the 
SNe at t < -1-25 days from by fitting Gaussian pro¬ 
files. From the figure we conclude that SNe 2011fe, 2012bm, 
2012cg and 2012cp show similar velocity evolutions, and the 
single measurement we have of SN 2012cu is also consistent 
with SN2011fe. The SNe 2012bl, 2012et and 2014J, on the 
other hand, are classified as HV. Note however that the 
velocity range for the full sample is still small enough for 
the scatter induced by the optical colour-velocity relation to 
land within the adopted uncertainty of 0.1 mag for B — V 


(Mandel, Foley & Kirshner 2014). Mandel, Foley & Kirsh 


(2014) further find no evidence for such a relation for 
the V — R and V — i. 

SN2012et is both classified as BL (although the pW& 
and pW7 values put its close to the border of the BL and 
CN groups) and HV, and such objects have been shown to 



-5 0 5 10 

Phase (days) 


Figure 12. Measured velocities of Si II A6355 for different phases 
of the SNe studied in this work together with reference objects. 
For SNe 2011fe and 2014J the velocities were calculated based 
on the spectra published by [Pereira et ah] | |2013[ ) , [Goobar et al.| 
) |2014| l and [Marion et al.j | [2015 . The grey region is marking the 
average SN la velocity range Folatelli et al.[[2012| [201^ . The 
dashed line marks the difference used to separate high and low 
velocity objects by [Wang et al.[|[2009[) at maximum. 


be redder than the bulk of normal SNe la (Folatelli et al. 


2013 Burns et al. 2014). However, the value of pWl for 


SN 2012et at —8 days is pWl= 82 ± LA, which is well below 
pVFl~ 150 A where the BL HV start to deviate in colour 
(see Figure 17 in Burns et al.|2014 i. 

M13 also determine that NUV-blue SNe show evidence 
of unburned C ll in their optical spectra while the remain¬ 
ing groups consist of objects both with and without un¬ 
burned C II detections. These results are also consistent with 
optical spectroscopic studies where SNe la with signatures 
of unburned carbon are bluer in their optical colours (e.g. 
Thomas et al.[ [2011[ [Folatelli et al.[[2012[ [Silverman fc FiL| 
ippenko 2012[ ). Unburned carbon can be detected through 
absorption attributed to the most prominent line in the op¬ 
tical, C II A6580, but the prospects of detecting this line 
decreases for epochs approaching maximum light (see e.g. 
Figure 11 in Folatelli et al.|[2012 ). The SNe 201 Ife, 2012cg 
and 2014J have all evidence for unburned material in their 
spectra (Nugent et al. 2011 Silverman et al.||2012 Marion 


et al.|2015 1, while the first spectra of the remaining objects 

were all obtained > — 7 days with respect to 

Finally, we can also compare the optical spectra of SNe 
2011fe, 2012cg and 2014J at the late-time nebular phase 
when the material in the SN ejecta is optically thin. In Fig¬ 
ure we present spectra of the three SNe at 260-280 days 
past maximum after they have been corrected for redden¬ 
ing. The spectra are dominated by forbidden lines from the 
decay chain ®®Ni —^®®Co —^®®Fe. The spectra of SNe 2012cg 
and 2011fe are remarkably similar while for SN2014J the 
[Co III] line at 6200 A is suppressed and the features at 
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Wavelength [A] 


Figure 13. Nebular spectra of SNe 2011fe, 2012cg and 2014J 
obtained with the William Herschel Telescope, the Keck and the 
Apache Point Observatory. The epochs of the spectra with re¬ 
spect to maximum is specified in the legend. The spectra of SNe 
2012cg and 2014J were corrected for reddening and smoothed us¬ 
ing a 5 pixel boxcar, and then scaled to match the spectrum of 
SN2011fe. The main features of the nebular spectra have been 
marked. 


7000-7400 A are redshifted with respect to the others. The 
two peaks identified as [Fe ii] A7155 and [Ni ii] A7378 (e.g. 
Maeda et al. ]|2010a|b| , are consistently blueshifted with re¬ 


spect to their rest frame by 800-1000 km/s in SNe 2011fe 
and 2012cg. In contrast, these lines appear redshifted in 
SN 2014J by a similar amount. This is in line with the find¬ 
ings of Maeda et al. (2010a I that supernovae showing a HV 
gradient ( Benetti et al.|2005 1, generally associated with HV 
supernovae ( Wang et al.|2009 1, demonstrate exclusively red¬ 
shifted nebular velocities, while LV supernovae preferentially 
demonstrate blue shifted nebular profiles. 

In conclusion, based on the optical spectroscopy avail¬ 
able both around maximum light and in the nebular phase, 
there are no observational signatures that would disfavour 
comparing the colours of the SNe to SN 2011fe for study¬ 
ing reddening together with the adopted uncertainties. In 
particular, while SNe 2011fe and 2012cg have slightly dif¬ 
ferent light curve shapes, they show remarkable similarities 
in spectroscopical features and evolution. However, the SNe 
2011by and 2011fe have also been found to be almost identi¬ 
cal in the optical while showing significant discrepancies in 
the UV ( Foley fc Kirshner|2013 Graham et al.|2015a l. Next 
we investigate the impact on the fitted reddening laws if we 
instead compare the measured colours to SN 2011by. 


6.2 SN 2011by as colour template 


SN2011by was discovered (Jin & Gao 20111 in NGC 3972 


and found to be a spectroscopically normal SN la (Silver- 
man, Ganeshalingam fc Filippenko|2013 1 with minimal red¬ 


dening ( Maguire et al.|2012p . HST/STIS UV observations 
were obtained at maximum (Program GO-12298; PI: Ellis) 
and using these data Foley fc Kirshner] (20131 and Graham] 


et al. (2015a I have shown that while SN 2011by is almost 


identical to SN2011fe in its optical properties it is intrin- 


Figure 14. The flux ratio between SNe 201 Iby and 201 Ife, 
and 2012cg and 2011fe at maximum after they have been cor¬ 
rected for Galactic extinction, corrected for reddening (in the 
case of SN 2012cg), and normalised between 5000-5500 A. We 
use the reddening parameters derived for maximum for SN 2012cg 
from Table and both spectra have been smoothed using box¬ 
car smoothening. These spectra are only shown for wavelengths 
where the signal-to-noise per spectral element exceeds 2.5. Some 
of the passbands have been plotted in the background to illus¬ 
trate which filters are affected by the difference in spectroscopic 
shape. The blue line shows Equation that is used to model 

the flux ratio between SNe 201 Iby and 201 Ife as described in 
Appendix [U] 


sically fainter in the UV. The two SNe are compared at 
maximum in Figure where the ratio between the spectra 
are plotted. Comparing this to the ratio between SNe 2012cg 
and 201 Ife at maximum, also shown in the figure, suggests 
that SN 2011by may be an even better reference SN than 
SN2011fe in the UV. 

Both SNe 2011by and 2012cg are fainter than SN 2011fe 
in the UV, which could explain why SN 2012cg is marginally 
redder in its UV colours than predicted by the fitted redden¬ 
ing laws. Two other SNe in the sample, 2012bm and 2012cp, 
show a similar behaviour for the bluest colour F225W — V. 

HST/STIS spectroscopy of SN 2011by was obtained at 
phases —10 and —1 days, which unfortunately does not al¬ 
low us to construct a SED covering the full range of our 
observations as we did for SN 2011fe. However, we can use 
the available spectrum at maximum and the Swift observa¬ 
tions of this SN to model the difference between SNe 2011by 
and 201 Ife in the UV. In Appendix [Cj we describe this pro¬ 
cedure and construct a model SED of SN 201 Iby under the 
assumption that it is identical to SN 2011fe at wavelengths 
A > 3000 A. 

Using this as the comparison SN, we can reht the red¬ 
dening laws for the SNe 2012cg, 2012cp and 2012bm for 
which we have reason to believe that SN 201 Iby may pro¬ 
vide a better comparison object. The htted F99 parameters 
at maximum are presented in Table The fitted parameters 
for SNe 2012cp and 2012bm change marginally compared to 
Table but we do note that reduced has decreased. How¬ 
ever, the most interesting result here is for SN 2012cg. Not 
only the htted reddening law parameters change, but also 
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SN A-range E{B — V) Ry tR R 


2012cg 


2012cp 


2012bm 


E{B- V) 
Rv 

xmax 

xV^ 


0 . 12 ( 0 . 02 ) 
o 7+1.4 

■i-' -1.0 
56080.5 (0.3) 
1.13(0.02) 
0.88 


0.34(0.03) 

q q + 0.5 
'^•^-0.4 

56081.1 (0.3) 
1.19(0.02) 
1.09 


0.45 (0.04) 
3 0’*"*^'^ 

'^•^- 0.4 

56018.3(1.1) 
1.24 (0.08) 
2.13 


Table 8. Best fitted parameters of the F99 law using the colour 
and SED model that has been adapted to match SN 2011by in 
the UV as described in §[^ The quoted uncertainties are the 
statistical Icr errors from the fit- 



Figure 15. The colour excess for SN 2012cg, similar to the upper 
left panel of Figure]^ with the exception that the measured and 
predicted colour excesses are based on the SED model described 
in §1^ which has been derived to mimic SN 201 Iby in the UV. The 
solid lines represent the best fitted F99 law with the parameters 
given in Table 


the time-evolution of E{uvm2 — V) becomes less significant, 
which can be seen in Fignre p^ This can be explained from 
Figure [nT| where the colour excesses of SN 2011by are shown 
with respect to SN2011fe, and a similar time-evolution for 
E{uvm2 — V) is shown. We note that the there is still a ten¬ 
sion between the bluest colours and the fitted reddening law 
which may suggest that the spectrum of SN 2012cg drops 
faster than SN 201 Iby at A < 2500 A. 

From the spectroscopic similarities between SNe 2011by 
and 2012cg as shown in Figure p^ and the similar UV-optical 
colours we find that SN 201 Iby appears to be a better com¬ 
parison SN for SN 2012cg than SN2011fe. This study also 
show how studying the UV-colour evolution of SNe la pro¬ 
vides a powerful tool for comparing predictions of CS dust 
scenarios with intrinsic colour variations. 


SN2012cgt 

UV-NIR 

0.11 

(0.02) 

3.8ll;5 

0.71 

OPT-NIR 

0.11 

(0.05) 

4.0l?J 

0.56 

SN2012et 

UV-NIR 

0.16 

(0.02) 

6 d 
+ 1 
GO 

1.06 


OPT-NIR 

0.15 

(0.09) 

<1.9 

1.72 

SN2012cpt 

UV-NIR 

0.34 

(0.03) 

q 1 +0.5 

0.83 

OPT-NIR 

0.36 

(0.04) 

q r» + 0.6 

1.03 

SN2012bmI 

UV-NIR 

0.45 

(0.03) 


1.68 


OPT-NIR 

0.41 

(0.05) 

q K+0.6 

0.70 

SN2012CU 

UV-NIR 

0.99 

(0.03) 

2.8+1 

0.47 


OPT-NIR 

1.00 

(0.05) 

‘ -0.2 

0.43 

SN2014J 

UV-NIR 

1.36 

(0.02) 

1 4+0-1 

-*-•^-0.1 

2.01 


OPT-NIR 

1.28 

(0.04) 

1 

2.29 


tSN colours were compared to SN 201 Iby 


Table 9. The best fitted F99 parameters using either only the 
optical-NIR colours or all colours (UV-NIR) while the and 
sg parameters were fixed to the values given in Tables and 
respectively. All colour data between phases —10 and +35 days 
from were used for the fits. 


6.3 The impact of the UV-data 

Another approach for testing the robustness of the reddening 
laws is to study how the fitted parameters are affected by the 
data at different wavelengths. We do this by comparing the 
fitted F99 parameter values for two different cases, where 
we use all measured colours (UV-NIR) and where we omit 
the UV data (OPT-NIR). A tension between the two cases 
could be a sign of either assuming an incorrect reddening 
law or colour reference SN. We keep and sb fixed for 
this study and use SN 201 Iby as reference for the SNe with 
red UV colours. The results are presented in Table 

The results are consistent for all SNe except for 
SN 2014J, where the addition of the UV data shifts the fit¬ 
ted values by ~ 2-3a. We fit a low value of Rv for this SN 
and, as shown in Figurej^ a power-law appears to provide a 
better fit in this case. In Table we present the results on 
the two colour sets when using a power-law relation for the 
two low-Rv SNe. While there is a slight, ~ la, tension for 
SN 2012et, the results are fully consistent between the two 
cases for SN 2014J. 

The major impact of adding the UV-colours is, how¬ 
ever, on the uncertainties of the fitted values, in particu¬ 
lar for low-reddening SNe. This is illustrated in Figure 
where the grey contour shows the 68 % confidence region for 
E{B— V) and Rv when the F99 law is fitted to the optical- 
NIR data of SN 2012cg, while the blue contour shows the 
uncertainties from the full UV-NIR fit. From the figure we 
conclude that extending the wavelength range significantly 
improves the constraints on Rv despite the high intrinsic 
dispersion adopted for the Swift and HST colours. 
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SN 

A-range 

Ay 

0 


SN2012et 

UV-NlR 

OPT-NIR 

0.30 (0.13) 
0.16(0.11) 

-1.7 (0.5) 
-3.3 (1.5) 

0.63 

0.82 

SN2014J 

UV-NIR 

OPT-NIR 

1.97 (0.10) 
1.93 (0.10) 

-2.0 (0.1) 
-2.1 (0.1) 

1.40 

1.02 


Table 10. The best fitted power-law parameters for the two SNe 
with low Ry values. The fits were carried out using either only 
the optical-NIR colours or all colours (UV-NIR) while the 4^“’' 
and SB parameters were fixed to the values given in Table All 
colour data between phases —10 and +35 days from were 

used for the fits. 



Figure 16. Contours showing the joint 68% uncertainty region 
for the fitted parameters Ry and E{B — V) for a F99 extinction 
laws for SN 2012cg in comparison to SN 2011by. The two contours 
show the fit results when the full data set (UV—NIR) is used (dark 
blue contour region) and when the colours based on HST and 
Swift UV observations are omitted (light grey contour region). 
All data between phases —10 and +35 days with respect to 
were used, and the best fit values for the two cases have been 
marked by black and blue dots respectively. 


6.4 Extinction from interstellar absorption 
features 

In §[^we discussed the constraints on the existence of cir- 
cumstellar dust that can be derived from non-detection of 
time-varying Na l D absorption lines. Na l D lines in high- 
resolution spectra are also commonly used as a proxy for 
reddening ( Poznanski et al!]|2011 l and the equivalent width 
of the lines have been shown to correlate well with extinction 
in the Milky May (see Munari fc Zwitt^|1997 Poznanski, 


Prochaska & Bloom 20121. However, Phillips et al. (20131 


found that of their sample of SNe la ~ 25 %, of which most 
have blue shifted Na l D features, have unusually high Na l 
column densities compared to the values expected from the 


Poznanski, Prochaska & Bloom (20121 relation. It is cur¬ 


rently unclear whether the high abundance of sodium is in¬ 
trinsic to the environment of some SN la or if the Milky 



Wavelength (A) 


Figure 17. Mg II at AA2796,2803 A and Mg I at AA2853 A 
marked by dotted lines in STIS spectra of four SNe with dif¬ 
ferent reddening. The continua of the unsmoothed spectra are 
normalised and offset for clarity. Above the respective spectra, 
the phases are indicated at which the spectra were taken. 


Way reference value is anomalous compared to other galax¬ 
ies. Other species, such as diffuse interstellar bands, have 
been shown to be useful for studying the ISM in SNe la host 
galaxies ( Sollerman et al.|2005 1 and in particular the diffuse 
interstellar band (DIB) feature at 5780 A may provide a 


better proxy for reddening of SNe la than Na l D (Phillips 


et al. 20131. In Table |11| we present equivalent width val¬ 
ues of Na I D and the DIB at 5780 A with corresponding 


reddening estimated by the empirical relations (Poznanski, 
Prochaska &: Bloom|2012 | and ( Phillips et al.|20]^ l, respec¬ 


tively, for the SNe of which we have obtained spectra or 
published data exists. Further empirical reddening relations 
exist for ISM features Ca ll H&K ( Murga et aI.|[20T5 l and 
Mg II ( Menard et al.|2008 l. These relations have so far not 
been extensively studied for their validity in estimating the 
reddening of SNe. 

Besides Na l D, narrow Ca ll H&K features, with respec¬ 
tive equivalent widths of 770 ± 30 and 579 ± 35 mA could 
only be measured from the first FIES epoch of SN 2012cg. 
These values are well above the range in which the relations 
in ( Murga et al.|[2015 l are valid. Phillips et al. (20131 also 
quote the detection of the DIB 5780 A which is consistent 
with the non-detection in the individual FIES epochs within 
the signal-to-noise. The most prominent absorption features 
in Na l D, shown in Figure [To| roughly agree with Vh from 
Table and are slightly blueshifted with respect to the ap¬ 
parent stellar velocity of ~ 448 km/s along the line-of-sight 
of SN 2012cg (see Cortes, Kenney & Hardy 20061. Additional 
features of the Na l D and Ca ll H&K profiles span a range 
from ~ 370-530 km/s, whereby a number of features are 
blended on the redshifted side of the profile. SN 2012cg is 
included in the analysis of Phillips et al. ( 2013| ) and is not 
counted in the group of SNe with anomalously high Na l 
column densities. However, the best fit E{B— V) in Table[^ 
for 2012cg highly disagrees with the value computed from 
the Poznanski, Prochaska & Bloom (20121 relation. 

Four epochs of high-resolution spectra of SN 2012cu 
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SN 

Na I D1 

Na I D2 

E{B- V) 

DIB 5780 

Av 

Mg II 2796 

E{B- V) 

E(B- V) 


(mA) 

(mA) 

(Na I D) 

(mA) 

(DIB) 

(A) 

(Mg II) 

(Phot.) 

2014J 

2558 (6) 

2831(9) 

- 

344 (4) 

1.8 (0.9) 

^ 4.6 (0.2) 

0.14(0.08) 

1.36(0.02) 

2012cg 

700 (9) 

1018(8) 

-*-•^-0.3 

85(5)1 

0.4 (0.2) 

^ 3.5 (0.2) 

0.08(0.05) 

0.11(0.02) 

2012CU 

849 (3) t 

925 (3)t 

. 7+0.4 
' -0.3 

- 

- 

- 

- 

0.99 (0.03) 

2012et 

- 

650(40)* 

Q O+0.3 

- 

- 

- 

- 

0.16(0.02) 

201 Ife 

27.4(1.1)* 

47.1(0.8)* 

i -0.003 

- 

- 

2.9(0.2) 

CO 

O O 

d d 
+ 1 

o 

o 

- 

201Iby 

- 

- 

- 

- 

- 

^ 2.3 (0.2) 

0.04(0.02) 

- 


and 


Patat et al. 


120131. 


Phillips et al. 

(2013 


Sternberg et al. 

(2014 


Maguire et al. 

(2013 


Table 11. Equivalent width values of Na I D, the DIB feature at AA5780 A and Mg II at AA2796 A for which reddening parameters can 
be inferred from [Poznanski, Prochaska &: Bloom] l (2012[ |, p^hillips et al.| | |2013[ l and |Menard et al.| l |2008l l, respectively. Upper limits are 
determined for Mg II measurements, since the doublet is blended in the STIS spectra. In the case of SN 2014J, the Na I D equivalent 
width extrapolates the relation by [Poznanski, Prochaska fc Bloom] | |2012[ l beyond any reasonable value for E{B — V). Best photometric 
E{B — V) values have been taken from Table |9| for comparison and the E{B — V) for 2011fe and 2011by is assumed to be negligible. 


have been described in Sternberg et al. (20141. No time- 


variability is detected in the deep Na l D absorption features, 
which likely consist of several blended components. The ve¬ 
locity of the Na l D lines of ~ 1130 km/s are in agreement 
with the H I velocity along th line-of-sight of SN 2012cu 
in NGC 4772 ( Haynes et al.|[?000| . We perform the photo¬ 
ionisation analysis described in Sectionj^on the Na l D dou¬ 
blet of 2012cu. The non-variation of Na l D1 indicates that 
there is no Na l between 6-10^® - 2-10^® cm. However, these 
constraints are weak since the absorption lines are close to 
saturation and not optically thin. The E{B — V) computed 
from the Poznanski, Prochaska & Bloom (20121 relation is 


~ 2a greater than the photometric value. Furthermore, this 
value must be considered a lower limit since the Na l D 
doublet of SN 2012cu appears to be close to saturation. 

SN2012et is included in the sample of [Maguire et ah] 


(20131 in which an Na l D2 equivalent width is quoted from 


a mid-resolution spectrum. The corresponding [Poznanski,] 
Prochaska & Bloom (2012 I relation E{B— V) is ~ Icr greater 
than expected. 


It is clear from Phillips et al. (20131 that Na l D do not 


serve as good reddening proxies for SNe la and that the scat¬ 


ter in the Poznanski, Prochaska & Bloom (20121 relation is 


underestimated. This appears to be the case for SNe with 
high reddening, such as 2014J, as well as low reddening such 
as 2012cg. Furthermore, many of the Na l D profiles do not 
appear to be optically thin, implying that the column densi¬ 
ties are rather underestimated, which would further increase 
the discrepancy between the |Poznanski, Prochaska fc Bloolii] 
(20121 relation and observations. There is thus motivation 


to investigate other ISM species, such as e.g. Mg ll, which 
are more difficult to observe and have not been studied in 
large samples so far. 


STIS spectra of SNe 2014J (from Foley et al. 20141, 
2012cg, 2011fe and 2011b3Q around the Mg ii doublet at 


AA2796, 2803 A are shown in Figure 17 The absorption fea¬ 


tures can be used to test the relation 


E{B- V) = (0.008 ± 0.001)IF^^ ®®=^°'^'^' 


^ STIS spectra of SNe 2011fe and 2011by were taken as part of 
Program GO-12298; PI: Ellis 


of Menard et al. (20081, where W is the equivalent width 
of Mg II at AA2796 A. The STIS spectra of SNe 2012bl, 
2012bm, 2012cp and 2012et are not included in the sam¬ 
ple due to the low signal-to-noise around the Mg ll doublet. 
Upon visual inspection of Figure Eli one can see that the 
Mg II lines are deep for all four SNe and the line ratios ob¬ 
tained from fitting a double Gaussian profile further indicate 
that the lines are not optically thin. In all STIS spectra, ex¬ 
cept those of 20life, the Mg ii doublet is strongly blended. 
Further, it cannot be excluded that the entire absorption 
is attributed to the Milky Way. We therefore measure the 
total equivalent width of the observed features and consider 
2/3 of the value to be the upper limit of the component at 
AA2796 A. The fraction is chosen from the line ratio of the 
Mg II doublet, if it is optically thin. The limits we obtain 
along with the corresponding E{B — V) computed from the 
relation above are shown in Table EH The Mg II doublet is 
resolved better in the 2011fe STIS spectra and the equiva¬ 
lent width can be determined from the Gaussian fits more 
accurately. The value presented in Table ED is the average 
equivalent width computed from all available 2011fe STIS 
spectra. In comparison to the photometric E{B — V), the 
[Menard et al.| ( [2008[ ) relation clearly underestimates the red¬ 
dening of 2014J, whereas no firm conclusions can be drawn 
for the other SNe due to the large errors and since blending 
with Milky Way features cannot be excluded. 


6.5 The origin of SN reddening 

While it can be expected that the observed SN la colour 
arise from a mixture of extinction by dust and intrin¬ 
sic SN la colour diversity, the latter is typically of order 
|A(i3 —1/)| < 0.1 (e.g. Nobili fc Goobar|2008 1. Thus, we can 
assume that the impact of dimming by dust will dominate 
with increased observed reddening and that the htted extinc¬ 
tion laws will be less suceptible to intrinsic SN la variations. 
It is therefore particularly interesting to compare the fits for 
SNe 2012cu and 2014J. They are the reddest objects in our 
sample, and at the same time, show very different total-to- 
selective extinction, Rv = 2.8 ± 0.1 and Rv = 1.4 ± 0.1 for 
SNe 2012cu and 2014J respectively. 

The evidence for material in the line-of-sight is strength- 
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ened by the fact that both SNe show deep Na l D absorption 
lines in their spectra (e.g. Sternberg et al.|2014 Goobar et al. 


2014 for 2012cu and 2014J respectively). Since CS dust is 


expected to give rise to steep extinction laws, such a sce¬ 
nario may offer an explanation for the low Rv of SN 2014J. 
However, no time-variability has been found for the Na l D 
absorption lines for any of the SNe we studied, although 


Graham et al. (2015b I report a detection of time-varying 


K I absorption for SN 2014J. When combined with the non¬ 
varying Na I D, these findings point to circumstellar dust 
at radii of rdust ^ 10^® cm due to the difference in ioni¬ 
sation energies. The detection of light echoes ( Crotts|[2015 1 
also points to the existence of dust at a comparable distance 
from the site of the explosion, 11 pc, i.e., rdust = 3-10^® cm. 

Further, we have studied the point spread function 
shape of SN 2014J in the HST A14 photometry to search 
for signs of dust at smaller radii but none could be detected, 
which is consistent with the lack of time-evolution of the 
SN 2014J colours with respect to SN 2011fe. For CS dust at 
rdust sj 10^® cm the photon time-delay exceeds the typi¬ 
cal time-scale of the SN la light curve and no colour time- 
evolution is expected as discussed in|Amanullah & Goobar| 

Other observations suggesting that the nature of the 
extinction of SN 2014J is primarily interstellar, is that the 
polarisation angle of the SN light is well aligned with the 
spiral structure of the host galaxy (Patat et al. 20151, the 
lack of thermal emission from the SN environment in near 
and mid-IR ( Johansson et al.|2014 l, and that the velocities 
of the observed multiple Na l D lines are consistent with 
measured Hi velocities along the line-of-sight in the host 
galaxy (Ritchey et al.||2015l. The latter is also the case for 


SNe 2012cu (Haynes et al. 20001. Even if dust does exist 
at rdust ~ 10^“ cm this is not expected to have any major 


impact on the measured Rv for such large distances (Wang 


2005). 


If the reddening of both SNe 2012cu and 2014J are dom¬ 
inated by interstellar extinction, the dust properties along 
the two llnes-of-sight must be significantly different, and the 
dust in M82 must be of a nature that has not been observed 
in the Milky Way (although see Brown et al. 2015 for 


a possible explanation). M82 is a starburst galaxy, and the 
dust properties have been studied by Hutton et al. (2014). 
They find that a steep Milky Way-like law is preferred over 
the Calzetti (20011 law, used for starburst galaxies and they 
derive a wavelength-dependence of the scattering within a 
central projected radius of < 3 kpc that suggests that ’’only 
small grains are entrained in the inner SNe-driven wind”, 
which is consistent with both the location (see Table and 
the characteristic of the observed reddening of SN 2014J. 

For SNe la, Rigault et al. (2013) found that when using 
Ha emission as a tracer for ongoing star formation, SNe with 
local Ha emission is on average redder and more homoge¬ 
neous, resulting in a lower brightness dispersion. These re¬ 


sults have been confirmed by Kelly et al. (2015 ) and Rigault 


et al. (2015) using UV data from the Galaxy Evolution Ex¬ 


plorer (GALEX). If the low Rv originates from small dust 
grains in regions undergoing intense star formation, we can 
expect to observe different local properties between SNe 
2012cu and 2014J. In Table fl^ we present the surface bright¬ 
ness of at the position of the SNe using an aperture with a 
2 kpc radius. From this we cannot draw any firm conclusion 


SN 

FUV 

NUV 

2014J 

26.5(0.2) 

24.9(0.1) 

2012bm 

26.7(6.1) 

25.9(2.4) 

2012cg 

25.1(0.4) 

24.4(0.2) 

2012cp 

24.7(1.3) 

24.6(0.7) 

2012CU 

26.6(2.3) 

26.0(1.0) 

2012et 

24.9(2.6) 

24.4(1.2) 

2012bl 

26.1(3.3) 

26.2(2.1) 


Table 12. Surface brightness (mag/arcsec) measured using a 
2 kpc radius centred on the SN position in the GALEX FUV 
and NUV filters. 


based on the local UV surface brightness of the host galax¬ 
ies. Even after we attempt to correct for host galaxy dust, 
following the recipe used in Rigault et al. (2015), we see no 
correlation with UV surface brightness and measured Rv- 
In Figure we show the fitted F99 parameters from 
Tablej^for all SNe together with similar measurements from 
the literature. Low values of Rv are not unique for SNe la, 
and has also been observed in the line-of-sight of other light 


sources such as gamma-ray burst afterglows (Fynbo et al. 
2014) and quasi-stellar objects (QSO, e.g. Fynbo et al.|2013 


Leighly et al.|20i^ . From the figure it is apparent that there 


are several observations with low Rv values for high extinc¬ 
tion SNe la, while SN 2012cu is both showing high extinction 
and a relatively high Rv - For low-extinction SNe it is by na¬ 
ture difficult to measure Rv for individual objects, but there 
are several statistical studies that obtain higher {Rv < 3) 
average values oi Rv compared to what is observed for high- 
extinction SNe la (e.g. Mandel, Narayan fc Kirshner||2011| 
Ghotard et al.pOll). The UV observations allow us to im¬ 


prove the sensitivity on Rv for individual low-reddened SNe, 
and it is interesting to note that we measure Rv discrepan¬ 
cies similar to what we saw for high-extinction objects. The 
colours of SNe 2012cg and 2012et are however a mixture 
of intrinsic variation as was shown by the comparison of 
SN 2012cg with SNe 2011by and 2011fe in particular. For 
these objects we therefore cannot draw any firm conclusions 
to what extent the measured Rv is solely a dust property. 
Note though that there are strong indications of dust along 
the line-of-sight for these SNe from Na l D absorption (as 
we have seen for SN 2012cg and for SN 2012et in |Maguire| 
et al.|[2M3l. Independent of the origin of the reddening we 


conclude that our sample of reddened SNe la with well ob¬ 
served colours, although limited, suggests that a range of 
Rv values can be expected for both low and high-reddened 
SNe la. 

Given the findings of M13 that all ”NUV-blue” SNe 
belong to the photospheric NV class while the ”NUV-red” 
consists of both NV and HV SNe and the colour-velocity 
relation for B — V ( Wang et al.|[200^ , we can also com¬ 
pare the derived extinction paramet ers wit h th e measured 
Si II A6355 A velocities. From Figures [l^ and 18 we conclude 
that both SNe for which we measure a low Rv (2014J and 
2012et) are classified as HV SNe while the remaining belong 
to the NV class. This is consistent with the statistical study 
of samples of SNe la by Wang et al. (2009), although our 
sample size does not allow us to draw hrm conclusions. How¬ 
ever, the fact that we present SNe with both high and low 
values of Rv in both high and low extinction environments 
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Figure 18. Fitted Ry and E{B — V) for the full phase range of 
— 10 to +35 days, where SN2011by has been used as reference 
for SNe 2012bm, 2012cg, and 2012cp and SN2011fe was used for 


Figure 19. The normalised colour excess E{X — V)/E{B — V) 
is shown using the same approach as in Figure]^ for the measure¬ 
ments of the three least reddened SNe. The errorbars are domi- 


the remaining SNe. Also plotted 

are measurements from the lit- 

nated by the adopted intrinsic colour uncertainties. Overplotted 

erature, where 9 SNe la from|Hernandez et al.|(|2000|l;|Krisciunas| 

is also the average F99 law (solid black) for the three SNe (fit- 

et aL| ||2006l |2007|l; |Benetti et al. 

l|2004||; |Elias-Rosa et al.|(|2006l 

ted by comparing the colours to SN2011fe for SN2012et and to 

2008|1; IHicken et al.| (|2009|| and 

Folatelli et aL| (|2010|l together 

SN2011by for the others) and the version of the SALT2 colour 


|et al.| | |2009[ | are included. The parameters measured by study- 
ing reddenin g of QSO by foreground galaxies me from[Elfasd6ttir| 


et al. (2006 I and Ostman, Goobar & Mortsell (2006 I 


law is only defined over the wavelength range 2000-8000 A. The 
grey regions show the RMS of the full sample with respect to the 
average reddening law for different wavelength bins. 


argues against that the origin of the discrepancy in observed 
reddening laws is solely due to intrinsic SN properties. 

One possible explanation to both of these findings could 
be that HV SNe primarily explode in environments with low 


Rv dust. Wang et al. (20131 and Pan et al. (20151 have 


shown that HV SNe are, on average, both located in more 
luminous hosts, and in brighter regions, closer to their host 
nuclei, compared the corresponding NV samples. Further, 


the results from Hutton et al. (20141 showed that a steep 


Milky-Way law is preferred for the central region of M82, 
and similarly the extinction towards the Galactic bulge also 
follows a steeper extinction curve {Rv ~ 2.5 |Nataf et al.| 
2013) compared to the Milky Way average. Both the dust 


and SN la properties are likely to be affected by the envi¬ 
ronment where they are produced and if HV SNe la orig¬ 
inate from younger and more metal-rich progenitors than 
NV SNe la ( |Wang et al.|[M^ , it is possible that the aver¬ 
age dust properties in the two environments are different as 
well. 

Studies of high-redshift SNe for cosmological applica¬ 
tions are unlikely to include highly reddened SNe like 2012cu 
and 2014J, but objects with similar reddening to SNe 2012cg 
and 2012et are common. In Figure we show the nor¬ 
malised colour excesses, E{X — V)/E{B — V), for the three 
least reddened SNe, where we have used the approach from 
Figure]^ to obtain the plotted wavelengths. Figurealso 
shows the average reddening F99 law and the SALT2 colour 
law ( |Guy et al. 20071 which has been empirically derived 
without attempting to discriminate between the different 
sources of SN reddening. Plotted here is the version used in 
the cosmology analysis from Betoule et al. (20141. Although 


the B-R bands, the SALT2 law starts to diverge from the 
data in the UV, wavelengths that becomes increasingly im¬ 
portant when SNe la at high-redshifts are studied. The lack 
of understanding SN reddening is one of the dominating as- 
trophysical systematic uncertainties when SNe la are used 
for cosmology (see e.g. Amanullah et al.|2010 Sullivan et al. 


2011 Suzuki et al.|2012 Betoule et al.|2014 | and further UV 
studies are needed in order to discriminate between different 
sources of reddening. 


7 SUMMARY AND CONCLUSIONS 

In this work we have presented and analysed broadband 
photometry of seven Type la supernovae covering the wave¬ 
length range 0.2-2 pm, obtained mainly from HST, Swift 
and the Nordic Optical Telescope. 

All SNe are found to be normal with respect to their 
light curves and optical spectroscopic properties. The SNe 
have colour excesses in the range E{B — V) = 0.0-1.4 mag 
and we have studied the reddening properties by comparing 
them to the normal pristine SNe la, 2011fe and 2011by. In 
particular, we have explored three different extinction laws: 


O’Donnell (1994), Fitzpatrick (1999), and Goobar (2008). 


the reddening laws are consistent at wavelengths covered by 


We have further studied the time-evolution of the fit¬ 
ted reddening laws and tested the consistency over a broad 
wavelength range. For the one SN, 2012cg, we do find indi¬ 
cations of time-evolution of the bluest colours and we have 
investigated if this could be explained by circumstellar dust. 
Further, we have also searched for time-varying absorption 
features, in particular Na l D, for SN 2012cg in multi-epoch 
high-resolution spectroscopy. The high-resolution spectra 
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were also used to estimate the extinction along the line- 
of-sight using the absorption of Na l D and Mg ll as dust 
proxies. The validity of using SN 201 Ife as a colour template 
was discussed and compared to the resulting reddening pa¬ 
rameters when using SN 2011by. 

We summarise the main findings as follows: 


(i) We observe a diversity in extinction properties charac¬ 
terised by a total-to-selective extinction, Rv, ranging from 
Rv ~ 1.4 to Rv ~ 3. 

(ii) The diversity between the fitted reddening laws for 
different SNe is significantly greater than the discrepancy be¬ 
tween the laws for a single SN. Although, we do find that the 
empirical F99 law provides a better fit for SNe with Rv ~ S 
while a power-law is preferred for the \ow-Rv SN 2014J. 

(iii) For the two reddest SNe, both with E{B — V) > 
1 mag, we observe significantly different values of Rv with 
Rv = 2.8±0.1 and Rv = 1.4±0.1, for SN 2012cu and 2014J 
respectively. Based on their Na l D profiles being consistent 
with the observed Hi velocity fields in their host galaxies, 
the lack of thermal emission and the observed polarization 
being aligned with the host galaxy plane for SN2014J, we 
conclude that the reddening is likely dominated by extinc¬ 
tion in the interstellar medium of the host galaxy. Further, 
pre-existing circumstellar dust at radii rduat < 10^® cm is 
expected to give rise to evolution of the colour excesses, 
E{X — V), but this is not observed for either SN 2012cu nor 
SN 2014J. 


(iv) We compare nebular spectra of SNe 2011fe, 2012cg 
and 2014J and conclude that all three are very similar at 
epochs -1-280 days from maximum, with SNe 201 Ife and 
2012cg being remarkably similar. 

(v) We compare the STIS spectra at maximum of 
SNe 2011by, 2011fe and 2012cg and conclude that while 
SN 2012cg is well matched to both SNe at A > 3800 A it 
is better matched to SN 2011by at shorter wavelengths. 

(vi) For SN2012cg we observe time-evolution in 
E{uvm2 — V) with respect to SN 2011fe, but when 
the object is compared to SN2011by the evolution is 
less significant. Circumstellar dust at rdust ~ 10^^ cm is 
consistent with E{uvm2 — V) but disfavoured by the lack 
of evolution of the redder colours. However, this example 
illustrates that observations over a wide wavelength range 
provides a powerful tool for separating the origin of colour 
excess evolution of SNe la. 

(vii) For low-reddening SNe such as SN 2012cg, adding 
UV data to optical and near-IR can decrease the statistical 
uncertainty of Rv by > 50 %. 

(viii) The two SNe with low Rv are found to have the 
highest Si ll velocities at maximum which is consistent with 
previous studies ( Wang et al.|2009 1. 

(ix) By measuring the equivalent width from absorption 
of Na I, we confirm the results of Phillips et al. (20131, and 
find that the extinction derived for SNe 2012cg, 2012cu and 
2012et using the Poznanski, Prochaska & Bloom (20121 re¬ 
lation is higher than what we observe directly for these ob¬ 
jects. 


Future studies, including a larger sample of low redden¬ 
ing objects covering a wide range in host galaxy morpholo¬ 
gies, will be necessary to robustly estimate the impact of the 
diversity in the wavelength dependence of reddening on the 
estimates of cosmological parameters from SNe la. The cur¬ 


rent work, clearly showing that there is significant range in 
the extinction parameters, and possibly also extinction laws, 
should motivate further multi-wavelength studies of SNe la, 
also including the UV. 
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APPENDIX A: SUMMARY OF SNe STUDIED 
IN THIS WORK 


The seven SNe la that we study in this work are briefly 
summarised below: 


• SN2012bl was discovered by the Chilean Automatic 
Supernova Search (CHASE) 42'.'9 east and 14'.'6 north of the 
center of the galaxy ESO 234-19 on Mar. 26.38 UT from an 
unfiltered image (Pignata et al.|[2012l. A spectrum (Prieto 


20121 was obtained with the 2.5 m du Pont telescope at the 


Las Campanas Observatory on Mar. 27.28 and object was 
classified as a SN la using SNID ( Blondin fc Tonry|2007 l. 

• SN 2012bm, located 10”95 west and 10.4 north of 
the centre of the galaxy UGC 8189, was discovered on 
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Mar. 27th by Puckett et al. (20121. A spectrum was ob- 


2008). 


tained on Mar. 28.05 UT and it was classified (Cappellaro 
et al.|20T2) as a SN la using GELATO (Harutyunyan et al. 


SN 2012cg was discovered (Kandrashoff et al. 20121 by 


the Lick Observatory Supernova Search (LOSS, Filippenko 
|et al.|[200ip on May 17.220 UT. The redder than normal 


SN was located 17”3 east and l'.'5 south ( Cenko et al.|20i2 l 
of the center of the Virgo Cluster member NGC 4424 in 
a region with many blue stars and disturbed dust lanes. 


Graur & Maoz (20121 used pre-explosion HST images of this 


well-studied galaxy to rule-out most supergiants as possible 
binary companions for the progenitor scenario. 

• SN2012cp was discovered 6” 7 east and l'.'2 south of 
the center of UGC 8713 on May 23.2 UT by the Puckett 


Observatory Supernova Search (Cox et al. 


jevic (2012) and Zhang, Wu & Wang (2012) obtained spec- 


20121. Milisavl- 


troscopy of the object on May 25.3 and May 25.7, respec¬ 
tively, and reported that it was consistent with a SN la be¬ 
fore maximum. 

SN 2012cu was discovered on June 14.6 UT by |ItagaM| 


respectively. The light curves are then used to calculate the 
(A — U)o colours shown in Figure 

As explained in §|^we also need a template of the spec¬ 
tral energy density of the unreddened object in order to be 
able to calculate the expected extinction for each passband 
for a given extinction model. In A14 the extinction for a 
given set of reddening law parameters was first calculated 
for all available spectra of SN 2011fe and the extinction for 
any given epoch was then obtained using spline interpola¬ 
tions. Similar to the light curve spline model, this approach 
is adequate for the UV, but due to the sparse spectral cov¬ 
erage, it will cause the model to deviate from the light curve 
for the redder optical bands and NIR at > 25 days past max¬ 
imum. In A14 this did not affect the final results due to the 
high extinction of SN 2014J and the conservative intrinsic 
colour dispersions adopted in the analysis. 

In this work we have developed a more accurate SED 
template by first resampling the spectral time series used 
in A14 to a resolution of 10 A, which are shown as the red 
spectra in Figure [Blj These were then used to calculate the 
missing wavelength and temporal elements through linear 


et al. (2012) 3'.'l east and 27'.'l south of the nucleus of the interpolation (shown in grey in Figure Bll, and finally the 


galaxy NGG 4772 and later classified as a SN la by |Marion| 


et al. (2012). Pre-explosion G/iondra observations are avail¬ 


able of the host galaxy and were used by [Nielsen, Voss fc] 


Nelemans (2013) to constrain the X-ray emission from the 


progenitor scenario. 

• SN 2012et, located 5”3 east and 0”8 north of the center 
of MGG +04-55-47, was discovered by [Rich et al.| ( p012[ ) on 
Sep. 12.057 UT, and later classified on Sept. 13.80 UT as a 
high-velocity SNe la (Dennefeld et al.|2012). 

• SN2014J was discovered by [Fossey et al.| (|2014|) in 


the nearby galaxy M82 on January 21.805 UT, and later 


classified by Gao, Kasliwal & McKay (2014). The SN was 
located 54" 


west and 21" 
the host galaxy. 


south of the ill-defined nucleus of 


APPENDIX B: A SPECTRAL AND COLOUR 
MODEL FROM SN2011fe 

SN2011fe is a normal SN la with excellent temporal and 
wavelength coverage with negligible extinction along the 
line-of-sight and is therefore a suitable object for reddening 
studies. In A14 we used the available UV-NIR spectroscopic 
time-series ( Pereira et al.|20i3 Mazzali et al. 12014 ) and the 
NIR light curves ( Matheson et al.|2012 ) to derive the extinc¬ 
tion law of SN 2014J. In this work we extend this template 
with UV photometry ( Brown et al.|[2M2 ) and optical pho¬ 
tometry at 25 to 40 days past B-band maximum. 

The light curves in the UV F225W, F275W, F336W and 
the optical B VRi bands are first constructed by synthesizing 


full matrix is mangled to match the light curves described 
above. The mangling was carried out by first calculating the 
ratio between the light curve fluxes and the corresponding 
synthetic SED fluxes. A function was then setup for each 
phase where the calculated ratios were adopted for the ef¬ 
fective wavelengths of the filters. For all wavelengths shorter 
(longer) than the bluest (reddest) point, the function was as¬ 
sumed to take the same value as in the extreme points, while 
we interpolate the function value for all intermediate wave¬ 
lengths, using the calculated ratios. The SED for the given 
phase was then multiplied by the function. The difference 
between the synthetic colours calculated from the mangled 
SED and the spline colour model was always found to be 
within < 0.05 mag for all epochs, and often even smaller. 
In Figure |B2| these differences are shown for six different 
colours. Since these differences are within the adopted in¬ 
trinsic colour uncertainties they will not have any significant 
impact on the fitted reddening laws. 


APPENDIX C: A MODEL FOR THE UV OF 
SN2011by 

The SNe 2011by and 20life have been found to be remark¬ 
ably similar in the optical, both around maximum and at 
nebular phases, while showing significant discrepancies in 
the UV ([Foley fc Kirshner 2013 Graham et al.p015a ). In 


the spectra from Mazzali et al. (2014) and Pereira et al. 


12013). A light curve for each filter is then obtained by fitting 
smoothed splines using SNooPy ( Burns et al.|20lT ). For the 
redder optical R and i bands, the sparse temporal coverage 
after 20 days past B-maximum does not allow an accurate 
spline fit of the second bump. We therefore also use the 


Figure pl4| the ratio between the spectra are shown at max¬ 
imum to illustrate this. Unfortunately, SN 2011by does not 
have the wavelength and temporal coverage that SN 201 Ife 
has, but we can create a crude template for SN 2011by based 
on the assumption that it is indeed identical to SN 201 Ife 
with the exception of the UV wavelengths. 

In Figure[^the, solid, green line is showing the function 


fix') = c 


data from Munari et al. (2013) to constrain the spline fit 


{x'* + 1 )^''"^ 


- 1 


+ 1 with x' = ——— (Cl) 


at > 25 days past B-max. For the Swift/ UVOT uvm2 filter 
and the NIR JHKs we fit spline models to the measured light 


curves from Brown et al. (2012) and Matheson et al. (2012) 


where x is the wavelength given in microns. In the figure 
the function is plotted for the parameters xo = 0.2495, a = 
0.0360 and c = 0.3359. 
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1 2 3 4 5 

Wavelength ) 


Figure Bl. The spectral template constructed from measured 
spectra of SN 2011fe. Sections based on measured data are shown 
in red, while the interpolated sections, warped to the light curves 
are shown in grey. The filters used for the warping are shown at 
the bottom of the figure. 


SN 201 Iby was observed by Swift in the UV (Milne et al. 


2013 1 and in the optical (Silverman, Ganeshalingam & Filip- 
penko|2013 1, and the colour excesses of these measurements 


with respect to SN 2011fe are shown in Figure [CT| Here we 
note that E{uvm2 — V) is showing a similar time-evolution 
as was seen for SN 2012cg in the upper-left panel of Figure]^ 
Under the assumption that SNe 2011by and 2011fe are 
twins at wavelengths A > 3000 A we can create an ap¬ 
proximate template for SN 201 Iby by multiplying the SED 
of SN 2011fe, presented in Appendix [B] and shown in Fig¬ 
ure Bl by Equation (Cll. Further, if we fix the parame- 



-10 0 10 20 30 40-10 0 10 20 30 40 


Phase [days] 


Figure B2. Examples of the difference between synthetic colours 
calculated from the SED template, and the corresponding colours 
from the light curve spline templates. 


ters xq and a to the values given above, we can fit c for 
each epoch with uvm2 observations by forcing the synthet¬ 
ically calculated Auvm 2 , based on the SED to match the 
observed E{uvm2 — V) in Figure [cl| (again, assuming that 
y 20 iife _ y 20 iiby^ fitted values of c are shown in Fig¬ 
ure [C^ and the synthetic values of Auvm 2 that these values 
correspond to are been plotted as black dots in Figure [CT| 
Further, we parametrise c(p) with a Gaussian, 


c(p) = 


\/2/kc 


exp 


i(^). 


+ b , 


as shown by the solid line in the Figure |G2| with po — 
-0.939, a = 1.406, b = 0.211 and cr = 4.992. 

To summarise, we create a crude light curve and SED 
template of SN2011by by using the template of SN2011fe 
from §[B|and multiply it with Equation (G11, where the pa¬ 
rameters, including the time-dependent c, are given above. 


APPENDIX D: RED TAILS OF THE HST 
FILTERS 

The HST filter transmission used in this work, shown in 
the upper panel of Figure |D1| suffer from red tails. That 
is, part of the flux observed for an object in these filters 
originates from wavelengths that are redder than the range 
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Figure Cl. The difference in colour between SNe 2011by and 
2011fe, where the ground based photometry has been adopted 
from [Silverman, Ganeshalingam & Filippenko] ( |2013[ l. Observa¬ 
tions for which the colour was obtained from measured V-band 
data are shown with filled symbols, while values obtained from 
a V-band template are shown with open symbols. The colour 
difference between the SNe can also be calculated synthetically 
based on Equation under the assumption that this function 
describes the flux ratio of the two objects, for all epochs, after 
adjusting for the difference in distance. This is motivated by the 
observation that these SNe appear to be twins at all wavelengths 
except for the UV. The parameter c in the equation was fitted by 
matching the synthetic and measured uvm2 — V for each epoch, 
where the former is shown by the black points for the best fitted 
values. 


where most of the transmission lies. As mentioned in A14, 
the red tails are relatively small for the WFC3/UVIS fil¬ 
ters, e.g. the transmission of F225W is roughly one part in 
100 000 at 5000 A compared to the peak transmission at 
~ 2250 A. Despite this, the red tails will have implications 
when reddened SNe la are studied. For example, from spec¬ 
trophotometry of SN 2011fe at maximum (thin, blue line in 
the middle panel of Figure Dll, we conclude that ~ 0.5% 
of the flux (dashed blue line) in the F225W filter will orig¬ 
inate from wavelengths > 4000 A. However, extinction will 
suppress the UV-flux signihcantly and after applying a F99 
reddening law with E{B — V) = 1.4, Rv ~ 1.4, ~ 60% of 
the F225W flux (thick, blue line) comes from photons with 
A > 4000 A. 

Another way to illustrate the effect of the red tails is 
shown in the lower panel of Figure [Pl] where the solid lines 
show the effective median wavelengths for an observation 
of SN20life at maximum through the WFC filters with 
various amounts of reddening applied. The effective median 
wavelength is dehned such that 50 % of the light originates 
from redder and bluer wavelengths respectively. Further, 
the wavelength range defined by the points where 10 % and 
90 % or the light comes from bluer wavelengths (i.e. the re¬ 
gion that includes 80 % of the light) is shown for increasing 
reddening. From the hgure we conclude that, F275W and 
F336W are less plagued by red tails than F225W, where 
the latter, in fact, will be probing redder effective wave¬ 


Figure C2. The fitted values of parameter c in after the SED of 
SN 2011fe has been fitted together with Equation to match 
SN 201 Iby uvm2 photometry. The solid line is showing a Gaussian 
fit. 


lengths of the SN spectrum than the two former filters for 
E{B - V)>1 mag. 

The fact that we are comparing different effective wave¬ 
lengths when calculating the colour excess between a red¬ 
dened and unreddened source is true for all broad band ob¬ 
servations. However, fitting extinction properties based on 
filters with red tails that span a wavelength range where 
the SED drops rapidly, also puts demands on knowing the 
passband throughput accurately even when the transmis¬ 
sion drops several orders of magnitude as for the filters in 
Figure m 

The expected WFC3/UVIS passband red tails were 
tested (Brown 2011| ) in the final round of thermal vac¬ 
uum tests of WFC3 prior to launch. They studied the fil¬ 
ter throughputs in four wavelength bins centered on 450, 
600, 750 and 900 nm and found deviations of up to ~ 40% 
between the measured and expected transmission. If the de¬ 
viations are taken as an estimate of the throughput uncer¬ 
tainty at these wavelengths, we can evaluate its impact on 
the extinction analysis. 

We carried out this test by first applying a F99 red¬ 
dening law with E{B — V) = 1.4, Rv = 1.4 to the SED 
of SN 2011fe described in Section |4.1[ The impact on the 
calculated extinction for the two bands F225W and F275W 
(the test of F336W was not carried out due to limitations 
in the setup as described in Brown|2011 1 can then be calcu¬ 
lated. We do this by comparing synthetic magnitudes from 
the original SYNPHOT throughputs, to magnitudes obtained 
after the filters have been multiplied by a spline function 
derived from the deviations for the four wavelength bins 
mentioned above. The impact on the calculated extinction 
was found to always be < 0.05 mag. Further, the worst case 
scenario can be estimated by letting all the measured devi¬ 
ations from Brown (20111 maximise the amount of red light 
entering the filter. The result on the extinction from this 
was <0.10 mag for t < 20 days with respect to B max and 
0.15-0.20 mag for t = 20-40 days, which is still within the 
dispersion that we have adopted for the intrinsic UV colours. 
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Figure Dl. The upper panel shows the filter transmissions, 
Tx(A), of the WFC3/UVIS passbands, multiplied with an ar¬ 
bitrary normalisation, C. The thin, blue, line in the middle panel 
shows the spectrum of SN 201 Ife at B band maximum, while the 
dashed, blue, line shows the same spectrum after it has been mul¬ 
tiplied by the F225W throughput. The thick, blue, line shows the 
same spectrum after it has been reddened with a F99 extinction 
law with E(B — V) = 1.4 mag and Ry = 1.4, similar to what we 
observe for SN 2014J. All three spectra have arbitrary normalisa¬ 
tion, C. In the plot we have also marked the wavelength region 
from which 80 % of the light originates from, defined such that 
10 % of light comes from bluer wavelengths, and 10 % comes from 
redder wavelengths. The lower panel shows the median wave¬ 
length, i.e. where 50 % of the light originates from redder and 
bluer wavelengths, for three passbands for SN 201 Ife at B max¬ 
imum, reddened with a F99 extinction law (Ry = 1.4) with in¬ 
creasing extinction. The 80 % region is shown in the same manner 
as for the middle panel. 


It is also possible to investigate the accuracy of the 
full HST/WFC3/UVIS SYNPHOT throughput by compar¬ 
ing spectrophotometry of objects, with known spectra, to 
WFC3/UVIS photometry. This is in fact the most direct 
test we can carry out of our method for fitting extinction 
laws since this is using the very same procedure. In figure [D^ 
we have plotted stars from CALSPEc]^ with the correspond¬ 
ing UVIS-photometry using the same analysis path as in 

® http://www.stsci.edu/hst/observatory/crds/calspec.html 
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Figure D2. Comparison between spectrophotometry of CALSPEC 
stars using the WFC3/UVIS filter throughputs, with the corre¬ 
sponding WFC3/UVIS photometry against the synthetically cal¬ 
culated B —V colour. The star P330E was used for calibrating 
WFC3/UVIS. 


Section [2.1 1 with a 0.4” radius. The agreement is in general 
more than sufficient for the purpose of this analysis, but it 
is unfortunate that data are not available for redder stars 
for which the impact of the red tail can be expected to be 
more significant. The star P330E should be discarded from 
the comparison since this star has been used for calibrating 
the WFC3/UVIS SYStem. 

However, we can carry out the same test for the red¬ 
dened SNe la that have been observed with both STIS and 
WFC3/UVIS. In Table Dl we present spectrophotometry 
at maximum for SNe 2012cg and 2014J, using the trans¬ 
missions of the WFC3/UVIS filters, compared to the corre¬ 
sponding WFC3/UVIS photometxj. For 2014J the spectrum 
and photometry were obtained on consecutive days, and the 
magnitudes based on the former were correct to correspond 
to the WFC3/UVIS epoch using the SN2011fe template. 
These corrections are <0.1 mag. 

For these two SNe, the discrepancies between two in¬ 
struments are within 0.1 mag, including F225W. This is the 
filter with the most significant red tails, and the photon 
wavelength distribution for F225W is significantly different 
for two SNe. The consistency between spectrophotometry 
and WFC3/UVIS photometry for two SNe that define the 
end-points of the reddening range in our study is showing 
that any errors of the filter transmissions will not affect the 
conclusions drawn in the paper. 

Further, the fact that we also observe a discrepancy of 
~ 0.1 mag for the F336W filter, which has minor red tails. 
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SN 

Filter 

STIS 

WFC 

bO 

F225W 

15.89 

15.98(0.05) 

CM 

F275W 

13.67 

13.61 (0.05) 

o 

CM 

F336W 

11.85 

11.74(0.06) 


F225W 

18.43 

18.50(0.04) 


F275W 

16.55 

16.55(0.04) 

CM 

F336W 

13.15 

13.02(0.04) 


Table Dl. Spectrophotometry of STIS spectra at maximum of 
SNe 2012cg and 2014J at maximum using the filter transmissions 
shown in Figure |D1| compared to WFC3/UVIS photometry of 
the same objects. For SN 2012cg the observations were obtained 
in consecutive HST orbits, while the data for SN 2014J was taken 
1 day apart. For this SN, we have calculated the magnitude dif¬ 
ference synthetically from the SN 201 Ife spectra described in Ap¬ 
pendix 1^ Since the colour evolution is slow between the two 
observations, any discrepancy in the filter transmissions will be 
of second order for this correction. 


suggests that the discrepancies we observe for F225W may 
not necessarily originate from uncertainties of these. 
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